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FOREWORD 


Two approaches for the modeling of turbulence exist. One is to start 
from an open hierarchy of profiles and higher order correlations. The other 
is to formulate a closed sequence of transport processes: evolution of mean 

profiles, eddy transport coefficients (viscosity, diffusivity, and damping 
or amplification rate), and relaxation or memory loss for the approach of the 
transport coefficients to equilibrium. The former approach uses an ambiguous 
closure, or an arbitrary hypothesis of the length scale of mixing. The 
latter approach is based on a theory of transport and the analysis of the 
spectral structure of turbulence. 

The author follows the latter approach and develops a new kinetic method. 
This intact report includes research reported quarterly during the first year 
of effort, September 1981 through August 1982, with a second report to follow, 
which will include research for September 1982 through August 1983. 
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CHAPTER A: QUARTER ONE 

Kinetic Basis of Cascade Transfer in Turbulence 

C. M. Tchen 

The Graduate Center and The City College of the City 
University of New York, N.Y. 10031 


ABSTRACT 

Among the transport functions which characterize the evolution 
of a turbulent spectrum, the cascade transfer is the only func- 
tion which describes the mode-coubling as the result of the 
nonlinear hydrodynamic state of turbulence. A kinetic theory 
combined with a scaling procedure is developed, to derive a 
kinetic equation of velocity distribution and the transition 
equations of path perturbations. These equations permit to 
formulate a transport theory of turbulence to investigate the 
eddy diffusivity, the eddy viscosity and the cascade transfer. 


Circumstances under which the transfer may appear in the form 



of a direct cascade or a reverse cascade are considered. 


I . INTRODUCTION 


The Navier-Stokes equation of motion, with a pressure 


gradient and a buoyancy as driving fields, can be used to de- 


scribe the micro-dynamical state of turbulance and to serve as 


the basis of a statistical theory. But the nonlinearities in 


the inertia and the driving fields will confront us from the 


start of the analysis. Various statistical methods have been 


proposed in the literature but were not too successful in sol- 


ing the closure of the hierarchy of equations. By changing the 


velocity function into a "particle" velocity as an independent 


random variable, the kinetic method is relieved of this diffi- 


A-2 


culty. What remains is a nonlinear interaction between fluid 


and "particle", with a diffusivity as an integral operator to 


represent the memory in the non-Markoff process. 


The kinetic theory of turbulence presents the following 


three problems: (a) the kinetic equation of turbulence, (b) 


the transport theory of turbulence for the determination of 


the eddy diffusivity, (c) the kinetic basis of the cascade 


transfer. In the derivation of a kinetic equation, it is im- 


portant not to fall into the Bogoliubov hierarchy. ^ ^ Its 

8 

closure has not been as successful in fluids as in plasmas. 


This diffuculty may have prevented further progress in the 


kinetic theory of fluid turbulence. We introduce a "scaling 


procedure" to avoid this hierarchy and the involvement of a 


system of equations of singlet and pair distribution functions. 


A- 3 



The eddy diffusivity which characterizes the kinetic equation 


is defined as the time integration of the Lagrangian corre- 


lation of field fluctuations. We develop a transport theory 


for the transformation of a Lagrangian correlation into an 


Enlerian correlation. This relationship, or the time-space 


transformation, has attracted many investigators by using the 


Q 

hypotheses of independence and normality. The correction 


has also been evaluated. The problem referred to a dif- 


fusion model, i.e. to a micro-dynamical state without a driving 


field. This model would exclude the more realistic aspects of 


turbulence, as presented by the Navier-Stokes turbulence, the 


shear turbulence and the turbulent motions in a statified me- 


dium. We feel it important to find a generalization that in- 


corporates the driving field. This generalization brings up 
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new difficulties. Our kinetic theory of turbulent transport 


helps in switching the driving field to a new role of advection 


in the phase space and in prescribing probabilities of tran- 


sition to replace the above hypotheses of independence and 


normality. 


In order to be qualified as a transport property, the 


diffusivity has to reach" a statistical equalibrium within a 


finite time, called the relaxation time. This finite time is 


obtained, when a memory-loss can be found. The scaling pro- 


cedure selects and organizes the necessary mechanism of mem- 


ory-loss by the turbulent dissipations. 


The mode-coupling governs the transfer of modes across 


the spectrum in both directions toward the high as well as 


the low wavenumbers. It appears in the form of the moment of 
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the fluid-particle interaction, in analogy with the nonlinear 


Landau damping in plasma turbulence. Formally, the fluid-par- 


ticle interaction can be expressed by a series of high order 


moments, as we know that a velocity distribution could be ex- 


panded in this way. It is nevertheless an insufficient way of rep 


resentation. This indicates that, with these limitations, a 


pure continuum method on its own may not be successful in ana- 


lyzing the mode-coupling, since it had failed in describing the 


Landau damping. 


II. MICRO-DYNAMICAL STATE OF TURBULENCE IN THE CONFIGURATION 
AND PHASE SPACES 

The hydrodynamical equations of the motion of an incom- 
pressible fluid are the equation of momentum 



(la) 


(V *-)i * i . \* . 

and the equation of continuity 

v- & = o y (lb) 

MV ** 

where the evolution of the total fluid velocity u (tyx) under 
the force E fcxj is governed by a differential operator 

L(t,*J s u v - (2) 

wv *«-> ' 

and V is the kinematic viscosity of the fluid. The force 

E = - jr v !° (3) 

- j ■** 

may be due to the gradient of pressure p in a fluid of density 
f t or may take a more complicate form involving the tempera- 


ture difference and the buoyancy force as in a strati- 


fied medium. In all cases, the force can be written in the dif- 


ferential form 


VEfcxj = VP: u{%x) u(t*y?= Z(t,x) (4a) 
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or in the integral form 


f(t*J ~ ^(t*/x ') J Ji(% *')j 


(4b) 


where 




(5a) 


is Poisson's integral operator, such that 


- - j £ i( %*’) . 


(5b) 


Both forms result from applying the condition of incompress- 


ibility (lb) upon (la) . Here and in the following the inte- 


gration is understood to extend from - &o to «>o . 


The hydrodynamical equations (la) and (lb) 


can be considered as describing the micro-dynamical state 


of turbulence in the configuration space t,x. This can be 


transformed into the phase space t,x,v, by introducing a 


S- 


function 


A- 8 


(6) 


N(t,X. s f[\/-U.(t,x)] 


which satisfies the equation of evolution 


Off- L)N(t,x,v) - 0, 


( 7 ) 


with a differential operator 


L (t,x,vj= v-p - yv 2 + E-'b , 


(8) 


in the phase space. The force can again be written in the inte- 


gral form: 


§(U) - h(K 


( 9 ) 


with Poisson's integral operator 




,/ i 


(* -X 'I 


and the source function 


(10a) 


'ifx' v') s VV : v V 

* { vv) % 


(10b) 


in the phase space. 


12 
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The micro-dynamical state in the phase space (7) possesses 

certain important advantages over that in the configuration 

space (la) in eliminating the nonlinearities. This is done 

by changing the role of the velocity as a random function 

into the role as an independent random variable. By the same 

token, the inhomogeneous equation of micro-evolution in the 

configuration space is transformed into a homogeneous equation 

in the phase space. We shall fully utilize these advantages 

in developing our statistical theory of turbulence. Of cours 2 , 

the micro-dynamics should still keep a nonlinear form in the 

phase space, as seen from E • d N in (7) . Fortunately the 

treatment of this nonlinearity can be delayed by considering 

to be provisionally known. After we have closed the transport 

equations, we can pick up the nonlinear relationship (4a) which 
then will be in its uncomplicate form of a consititutive relation 
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The equivalence between the two spaces (1) and (7) , can be 


easily verified by taking the moments 

I , J<A\/ v N = u. (t,x) . 

of (7) . It is then seen that the equation (7) in the phase 


space actually corresponds to the simultaneous system of dy 


namical equations (la) and (lb) in the configuration space. 


III. THE BBGK HIERARCHY OF TURBULENCE 


If \/ is a stochastic variable by the stochastic behavior 


of LK.(t,x) , the function (6) can be called a distribution 


function, denoted by 

retaining nevertheless its stochastic character, with a 


(ID 


( 12 ) 
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micro - dynamical state (7) . The function (12) and 


the evolution (7) contain fluctuations with all the minute 


details which are unnecessary in a statistical study. A 


coarse-graining procedure which eliminates all stochastic 


fluctuations is the ensemble average 


A * O 


(13) 


called "global average". A fluctuating function subject to 


this average becomes deterministic, e.g. 


A N(t, x, v) * f(t v) , 

W V« I M M 


(14) 


and the deviation from the average is the fluctuation 

i~f - f ' (i5a) 

or, in terms of the operator notations, 

I -A = A , (15b) 


where "1" is the "unit operator", A is the average operator. 
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and A is the fluctuation operator. 


Upon applying the average operator A to the equation 


(7) of micro - dynamics, we get 


(\ + L) *)--A e(t*J ■ } fex.vj 

- - ? 7 |V> 2 . (16a) 


or 


(\+Lj^ ft, xyj = -^ ft x/ <; ^ x>'J ^ > */' *,v)j 


(16b) 


when use is made of (9) and of the notation 

* ($(ty,sS‘) f(t, *,''}) • (17) 

The functions (14) and (17) are called the singlet and the 


doublet distributions, respectively. 


In (16b) we see that the evolution of the singlet dis 


tribution depends on the doublet distribution, and that the 


latter evolution will expectedly depends on the triplet 
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distribution 


x, ■y) h (f(t <«■) ■ (18) 


and the sequence continues to form the BBGK hierarchy. This 


hierarchy is well known in the kinetic theory of plasma tur- 


bulence, and is seen here to reappear in fluid turbulence 
with the difference that is defined by (10b) in an 

incompressible fluid, while A (x,v)~ I in plasmas. l”4/8 

The closure of the hierarchy constitutes a difficult 


problem. It has been attempted by finding a small parameter 


in plasmas, but such a parameter is lacking in fluids. This 


may explain why so little progress has been made in fluid 


turbulence by the kinetic approach. The analysis of the spec- 


tral structure of turbulence requires the Fourier transform' 


ation of a correlation function, and therefore the second 
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moment of the doublet distribution. Thus this kinetic problem 


calls for the determination of both distribution functions. This 


obviously would pose a tremendous task. We shall devise 


a kinetic method, in which we can avoid the BBGK hierarchy and 


the doublet distribution function. To this end, we shall resort 


to the scaling procedure of fluctuations, which we shall de- 


scribe in the following. 


IV. THE SCALING OF FLUCTUATIONS 


The BBGK hierarchy which begins with the singlet distri 


bution (16) does not give any statistical information unless 


the doublet distribution is known. This means that 


a minimum of two equations of the closed hierarchy will be 
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needed to describe the fluctuations statistically. In order 


to avoid the BBGK hierarchy completely, we shall scale the 


fluctuation 


a - a l °\ a 


(19) 


into a macro-group E°' and a micro-group E 


The scaling is performed by writing the fluctuation op- 


erator 


A = /4 < °’V A ' 


(20a) 


,(o ) . ^ 

into two components n , n , which select the macro- 


group and the micro-group, respectively, and are called 


"scaling operators". Since (15b) is known as the "Reynolds 


decomposition" in turbulence, we may consider the scaling 


(20a) as a simple extension of this decomposition. 
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For the sake of convenience, we introduce 


A.-A + A' 


so that 


A'~ i -A. 


The two groups E ^ and E' have their durations of 


correlation and t' , with the inequalities 

c c 


tw ? T ' 
c c 


stating that the micro-group c is more random than the 


macro-group E'' . Their respective portions of the energy 

to- 

f( £ . ( ° >2 ) = f k dk' s(a') , j°Jk'S(k') i 

(0 h 

are clearly separated by the wavenumber variable tC i n the 
spectral distribution S(fe) of E - fluctuations, although 


the individual groups E and C may overlap in their 


Fourier compositions. Clearly, the two portions add to 
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(22b) 


t<£ 2 > = S(kj ■ 


Similarly, we scale the velocity fluctuation 


2 =u c °U u' 


Co) , 

into a macro-velocity U. and a micro-velocity U of 


energies 

i <^*>=pk'S(k') 

*0 k 

in the spectral distribution F(k.) of U - fluctuations. 


The two portions also add to 


/ z\ fk~co 

<“>- dk'F(k) 


The same scaling applies to the distribution function 


-V 

in the superposition 

'v / * /V 

f - f (> * i 


of a macro-distribution ^ and a micro-distribution -j' 


It is to be remarked that the scaling permits the 


(23) 


(24a) 


(24b) 


(25) 
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derivation of the spectral distributions by simply 


differenting (22a) and (24a) with respect to 71 , while this 


was n ot the case with the hon-scaled energies (22b) and (24b). 


Thus our statistical problem amounts to finding a kin- 


etic equation of the scaled singlet distribution only, instead 


of both the singlet and the doublet distributions in the un- 


sealed treatment of the BBGK hierarchy. 


It is to be stipulated that the decomposition into the 


scaled groups (20) differs from the decomposition into the 


Fourier components or a group of such components. The latter 


decomposition is a mathematical transformation without adding 


any physical concept to the treatment, while the former de- 


composition injects a scaling procedure of coarse-graining to 


A (®) 

the stochastic process by the operator n or the operator 
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• Consequently, the scaling procedure dis- 


tinguishes between a macro-process representative of the evol- 


ution of a macro-field and a micro-process representative of 


the transport properties as shaped by the more random micro- 


fluctuations in the medium in which the macro-field is propa- 


gating. Well defined probabilities of transition, that are 


self-consistent with the once given micro-dynamical state of 

A (p ) ~ / 

turbulence, will define the scaling operators A and n 


V. KINETIC EQUATION OF TURBULENCE 


The micro-dynamical state of turbulence is described by 


the dynamical state (7) in the phase space. We apply the scaling 


operators /4°^ and A ' and decompose (7) into the system 
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(26a) 


(\ + A'L)jYU,yJ = f°>) , 


(26b) 


which describes the evolution of the macro-distribution 


under the statistical effects of the more random micro-fluc- 


tuations E' ^ / while the ^ fluctuations are ex- 


cited by the gradients of the mean and the macro-distributions. 


according to (26 ). All operators including the scaling oper- 


a(°J a ' 

ators n , n and the differential operator L refer to all 


functions which follow them. 


If we disregard the group-coupling between 0) and f " 
and merely express E in terms of according to (9) by 


A , we would again fall into the BBGK hierarchy in the form: 


~ ^ 


with 


(\ f L J V J - */*', s'jj *(*,' y' ’jfCt *y')f '(t, x, *)J 

3 -5 - A; yy') I a6 <; *j f°{ t * / ' V *. yjJ 

f it x'.v'j xy) 2 A (o) j '(tx ; A 


A- 21 



by definition of the scaled doublet distribution function. 

In order to avoid this hierarchy which is equally undesir- 
able, we call on the group-coupling between the two equations 
(26a) and (26b) . For this purpose, we rewrite (26b) in the form 

(VL)f' = - At < 27 > 

by the use of (20c) . We integrate to find 

f = A 

* - f dx U(t y t-zJE (t-x) • -h ^(t-z)] • (28) 

J 0 

The left hand side of (27) represents the total time derivative 

Tt E 1 + L <29) 

along the trajectory of the evolution of (t/X,\/) in the 
phase space, and the right hand side is a driving field. The 

solution (28) is the time integration of the driving field in 

A-22 
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the Lagrangian representation denoted by 


t 1 = Bit- % #«j] -n- i z xJt rj Yfi _ z) • ,30, 

Here x(t-T) is the position along the trajectory at the 


time t -~C , when it is known that at time t, the position is 


x and the velocity is v . The trajectory is governed by: 


dx (t 1 ) 


dt. 


dV(t, ) 


V(t,) , — 

- x dt. 


= E(t,) , 


Ola) 


where V(t,) is an intermediate variable in 
*•* “** 


d 2 x(t x ) (31b) 

— = E(t.) . 

d ti 2 

The dynamical equations (31) governing the trajectory are called 


the micro-dynamical equations of transition. The operator 


is called the propagator, or "Lagrangian operator". The last 


term of (27) does not contribute to the solution, by AA^O • 


The stress i n the phase space — is thus found to be 



( 32 ) 
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where 



( 33 ) 


is the eddy diffusivitv. Hence, upon substituting the stress 


into the equation (26a) of macro-evolution, we obtain: 

■ <34) 

In the evolution of the macro-distribution (34) of scale t , we have 


c, 

(35) 


and a fortiori 


t » T ' 

c 


(36) 


from 21) , so that ]D from (33) becomes a large-time diffusivity 

' nt 

-s 

00 


r>'» A (o) j Jr E'(t,*)XJ(t,t-T)ETt-T) y 

"Z V/i " ^ 


(37) 


that asymptotically will take a deterministic form: 


- I 0 ^ r ( e '(*' E '(t-rjy , 


(38) 
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on account of the small correlation time (36) and the deter- 


ministic character (22a) . Under this circumstance, we reduce 


(34) to the following kinetic equation 

It has the form of a Fokker-Planck equation, containing a 


memory as described by the diffusivity operator l . By 


the group-coupling and the scaling procedure, the kinetic 


equation (39) is explicit in -f 7 , avoiding the BBGK hierarchy. 


KINETIC THEORY OF TRANSITION 


The transport theory for the determination of the diffusi- 


vity (38) requires the examination of the Lagrangian correlation 

<£rw t ? *&-*)]) 


( 40 ) 


A-25 


The trajectory in the phase space 


>(t-r) = Y - vz f 2(-t) 


( 41 ) 


contains a fluctuating path 


in 


in the retrograde tran- 


sition, i.e. in a time interval — X , and is self-consistent 


with the same micro-dynamical state (31) as that of the kinetic 


equation of velocity distribution (39). 


The fluctuating path %(~t) is a random function, for 


which we will develop a kinetic theory of transition. 


The first of the micro-dynamical equation of transition 


(31) can be written in the form 


]*■*!> -°- 


u 


(42) 


called the "master equation of transition", where the velocity 


Vjf-t ) has been defined by (31) , and 

f(- T,l) = *[i -{(*)} 


(43) 
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is a probability of transition in the microscopic description. 
It is a function of the random variable X • an< ^ has the 


moments : 


Ul P= /, Uf. l P = M- v ) ■ 


(44) 


It is not difficult to verify, by means of the moments (44) , 


that the master equation (42) will indeed reproduce the micro- 


dynamical equations (31) . 

We apply the scalings A , A 0 j A to the master 


equation (42) and derive the following kinetic equations of 


transition 



and the probability of transition in the micro-scale 






(45) 


(46) 


( 47 ) 
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where U is a drift velocity. The details of derivation. 


which are similar to those used in the derivation of the kin- 


etic eauation of velocity, distribution, will not be repeated 


here. 


The coefficient of diffusion, as defined bv 


is ' ^ 


is calculated from the micro-dynamics (31) of the path fluc- 


tuation, written in the new form : 

X 

}(-x) Vtf) . . 

- o - 'o 

We find the dispersion 

(t(-x)M- *h '(fr-f TVt 


and the diffusion coefficient 


=r T E(t-T 'jy 


By scalings, we have: 
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(51a) 


4‘tvfW) * t 

K'/t-c) - £T 3 (eHt) E ( ° J (xj) 

•V 

(rr-^rr-t^ s j-t* i>' 

« T 2 !)' . . (51b) 

The kinetic equation of transition, in the non-Markoff 
process of scales (45) and (46) , is thus derived from the 
master equation (42) in self-consistency with the originally 
stated micro-dynamics of turbulence, in the form (1) , (31) , 
or (49) , and therefore should not be confounded with the com- 
monly known transition equation which was derived phenomen- 

13 14 

ologically in the hypotheses of a Markoff process. ' Note 
also that the kinetic equation of transition determines a 
retrograde transition of oath ^ 
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VI. TRANSPORT THEORY OF DIFFUSIVITY 


By the Fourier decomoosition, we transform the formula 


(38) of diffusivity into 


J> V 4xJ^k'(<r'[k', c 1 


’■ k'-vT 


(52) 


with the use of (40) and (41) . The details of the transform- 


ation are elementary and have been omitted. F or the sake 


of abbreviation, we have introduced the Fourier function 

V-'U'jM]* *Xk') c l -■ - ( - Vj 


(53a) 


and the fluctuating intensity 


* f<U' -O'(h') 

^ ^ V Ifc MU* 


(53b) 


of spectral tensor 


<-*?«> -s(r(‘’W» 


(53c) 


When the Fourier transform is truncated within a length inter- 


val l )( , the truncation factor is , where is 
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the number of dimensions. 


It is to be noted that the integral 


/ U 7 ' V M-rJ] \ 

with a perturbed path , constitutes a Lagrangian corre- 


(54a) 


* 

lation of the E - field fluctuations, while the integral 


Uj < vza : k-vix 


(54b) 


without path perturbations is the Euler ian correlation. Hence 


the distinction and the relation between the two correlations 


reside in the path perturbations in the Fourier function 


<sr‘> * U(k.') 
^ •» 


(55a) 


This function can be decomposed into two parts as 


n -) (a . + a y «, • c - 

w 

= <V r > * <<r' ) . 


(55b) 
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since A -f- A ' > with 


<y;>- 

■* 

(v~x>* (^i'(k’) 


(56a) 


(56b) 


For the analysis of the first part (56a) , we calculate: 

4 e - ik ' J -(- T) s ui 
„'e C W- X) •J 4 JL P'f-% 


*'e. *** 3 r 4 x 

A JjM ^hhA) - lAc , 


from the Fourier transformation of the kinetic equations of 


transition (45) - (47) , so that we find 


<<r'> = <V<iy> /■£ 


with the notations: 


H* l.-i -L' 

J> 14,0 


^ 5 **)>[- k'k'tjlt t' Kfc ')] 


ik'- u. V 


A-32 


■mil 1 ii ■ 



The second part (56b) contains the fluctuations already 


found in (57) , and can thus be written in the form: 


<^;> _ -ck- 

j ?* - x 


i 




(60) 


The degenerated function 


fc-v-f) - - ten 

J 


can be considered as a time-integrated flux, so that, by com- 


with 


bining (60) and (61) , we obtain 

(z^y- r^xk'jyH , 

? * > ' ■ ^ r-Tjj.'wy = r *-• ' * 'f 

in the transport of s' of a path and a diffusivity D 


(62) 

(63) 

/ 

s * 


The results (58) and (62) for the two parts of the Fourier 


function (55b) add to 


(tU: In) - ('- r) <f(k')> h . 


(64) 
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and yield a diffusivity which can be written in the form: 


x>' - Uk- QKvj 


( 65 ) 


This expresses that the diffusivity is endowed by the spectral 


tensor and approaches its equilibrium at a relaxation 


time 




( 66 ) 


where H and have been defined by (59) and (63), re- 


spectively. 


Note that, since the relaxation time (66) contains 3) 


in the integrand, through £ and from (63) , (59) and 


(51b), the form (65) is the integral equation for the deter- 


mination of D' . 


We conclude that the Lagrangian correlation of the 


E - field fluctuations is expressed in terms of the spectral 
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function by the following formula: 


<e '(t*J V( t t-*J eJt-xpMk. > fc) a U'))M(x k;J , 


(67a) 


with 


M (c. k'y) Sr 


e L - - T [i - ?(k j] Hj-x, k'J, 


(67b) 


from (64) - (66). 


VII. KINETIC BASIS OF THE CASCADE TRANSFER 


We write the kinetic equation of velocity distribution 


(34) in the following explicit form: 

\ f A- zb* *7 tm) + *X- ? ^ 

- i-vj-if 0 *-*)} ■ 


( 68 ) 


We take the moment of (68) . The zeroth moment gives the 


equation of continuity 


V- = 0 


(69) 
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and the first moment gives the equation of momentum in the 


form: 


- (p) „ / - {*) (o) - \ 

\ u. + V. [ u.. ti ; +■ Lt: u- f-A u,- u . ) 

1 c r L j - / w 7 


(o) — 


with a collisionless interaction 


( (0) (?) (0) ) £-(°J 


C L 


l 


(•) 


( 70 ) 


(71a) 


// * fa}' 

t the production ) ' 


The advection 


jS ^ / /oj (fi) I — 2 (d 

the inertia n^licU-J , the viscous damping V V Lt ■ 

(/ ‘ y 0 

-(o) 

and the driving field c are all the regular terms recog- 

t 




nizable by isolating the macro-group from the hydrodynamical 


equation (la) which describes the micro-dynamical state of tur- 

T C°) 

bulence. The interaction J is new, and can be verified to 


T 

vanish if the memory should be absent. Since J can be 


identified as a scaled Reynolds stress from the micro-flux of 


momentum, i.e, 


f:’ - 


- 7 . [A } u. u' ) 

J [ ‘ J J ' 


(71b) 
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the memory is indeed essential to turbulence . 


From the mixing-length hypothesis/ the momentum flux of 


the gradient type, i.e. 


u: 


- k: ^ 

tx ~~i> X 
/ A. 


(72) 


has been established, giving an interaction in the forms 

J ( ‘> - 2 . K' . (73) 

‘ ’<1 n, \ 

The transport coefficient is called the eddy viscosity. 


The comparison between the above two forms of interaction 


suggests a means of isolating this eddy viscosity, provided 


the transport of the gradient type (72) is valid. On the con- 


trary, when the gradient transport is not valid, as is the 


case with certain stratified media, the kinetic interaction 


(70) could determine the correct type of transport appropriate 


to the problem. 
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Upon multiplying the momentum equation (70) by 


(P) 
IA . 

L 


and averaging, we obtain the energy equation: 


JL'b 
2 t 


= p ( ^ + Q*>- 7« - £ W - ti 


( 74 ) 


describing the rate of change of the kinetic energy, as governed 


by the following transport functions: 


production ^ . 


coupling 


cascade 

transfer 


> *v 


(75) 


viscous 

dissipation 


flux transport 




Among these transport functions, only the mode-coupling 


function, or the cascade transfer function 

7 

= “/wy w? 


(76) 
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possesses a structure that is characteristically kinetic, as 


involving the correlation 




( 77 ) 


(o) 

between the fluid velocity U and the singlet distribution 


f 


(*) 


. This coupling is analogous to the fluid-particle 


interaction which was considered to be the origin of the non- 


linear Landau damping in plasma turbulence. 


In (67). we have derived a formula for the conversion of a 


Lagrangian correlation into an Eulerian correlation (or its 


Fourier transform). . We apply the same rule for the Lagrangian 


correlation (77), and obtain: 


/p)(+ I #i T )\ fjiL" V / L" )\ u/ ) 


(78) 


where ft/j is defined by (67b) 


Now when this Lagrangian correlation is substituted into 
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(76) , we find the transfer function in the form: 




( 79 ) 


Noting that 


??} 1 


is an operator of integrations 


with respect to k and T , as prescribed by (65) and (66) , 


respectively, we can rewrite (79) in the explicit form, as 


follows: 

x v [ m(t, tv; vj> . <«°> 


VIII- CONCLUSIONS 


By means of a scaling procedure, we have developed a kinetic 


equation of turbulence and a transport theory 


of eddy diffusivity. The transport theory is based 


upon the transformation of a Lagrangian correlation into an 
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Eulerian correlation by means of transition probabilities 


which we have found from a kinetic theory of transition. 


The transfer function describes the mode-coupling, or the 


cascade transfer across the spectrum. It appears in the form 


of a fluid-particle interaction, and is therefore a kinetic 


phenomenon like the Landau damping. A memory should exist 


in this interaction and is generated by the Lagrangian corre- 


lation between the fluid velocity and the "particle" velocity 


distribution. However, a memory-loss should be present for a 


collisionless dissipation to yield a finite relaxation time. 


A preliminary examination of the expression of the cascade 


transfer reveals that the cascade process can persist in the 


direction of high wavenumbers , so that either a transport along 


the gratient, or a collisionless dissipation in the explicit 
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form of a Landau damping can appear. With this distinction, 


an eddy viscosity or a rate coefficient of Landau damping can 


be separately analyzed. In other circumstances, e.g. in the 


atmospheric turbulence with a stable stratification, a reverse 


cascade may appear, involving a transport counter to the gradi- 

15 

ent. The reversal is associated with the "gap” phenomenon. 


An analytical expression of the cascade transfer of this kind 


will clarify the gap phenomenon and the concept of "negative 


viscosity" . 
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CHAPTER B: QUARTER TWO 


Kinetic Theory of Turbulent Transfer with Double Memory-Loss 

C.M. Tchen 

The Graduate Center and The City College of 
The City University of New York, N.Y. 10031 


Abstract 

The transfer function governs the mode-coupling in 
strong turbulence. It is investigated here by means of a 
kinetic equation of turbulence, in consistency with the 
hydrodynamical system that describes the micro-dynamical 
state of turbulence. The equation of spectral balance is 
obtained by taking the moment of the scaled kinetic 
equation for the one-point distribution without the need 
of the two-point distribution. It is found that the trans- 
fer function is governed by two Lagrangian correlation 
functions and therefore by two memory-loss functions. 

The first correlation relates to two field fluctuations 
and finds the diffusivity as an integral operator. This 
perpetuates the memory to the second correlation which 
contains its own memory from the fluctuations of velocity 
and distribution of velocities. This correlation is known 
as the "wave-particle interaction" in plasmas. The com- 
petition between the two memory-loss functions and the 
cut-off of the memory form the essential basis for the 
closure of turbulence. We find two forms of transfer 
functions, for small and large scales, respectively. The 
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first form belongs to the gradient transport and isolates 
an eddy viscosity. The second form belongs to a collision- 
less damping without gradient, and isolates a rate coef- 
ficient of collisionless damping. The transfer function 
for the large wavenumbers is applied to the derivation of 
the Kolmogoroff law of turbulence with the analytical de- 
termination of the numerical coefficient. 

1. INTRODUCTION 

The micro-dynamical state of fluid turbulence is 
usually described by a hydrodynamical system, for example 
the Navier-Stokes equation of motion. This equation is 
nonlinear in the velocity function and is inhomogeneous 
by having a driving field which may be the pressure 
gradient. On the other hand, in the kinetic representation 
the velocity is an independent variable and does not cause 
a nonlinearity, while the driving field forms a part of 
the advection and does not warrant an inhomogeneous 
equation. It is known that a homogeneous equation for the 
micro-dynamical state of turbulence can more explicitly 
give an analytical specification of the path perturbations, 
since the homogeneous equation contains a detailed differ- 
ential operator of perturbations and therefore will yield 
by inversion an exact evolution operator. 

Unlike the kinetic theory of gases, the kinetic theory 
of the spectral structure of turbulence usually requires a 
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two-point distribution function, since the spectral func- 
tion is the Fourier transform of the two-point correla- 
tion function of velocity fluctuations. We see that such a 

procedure would involve a hierarchy in the form of the 

1-7 

Bogoliubov hierarchy . In order to confine ourselves to 
a single kinetic equation of one-point distribution and 
thus avoid the Bogoliubov hierarchy, we introduce a pro- 
cedure of scaling into a mean group and groups of macro- 
and micro-fluctuations, representing three processes of 

transport: macro-evolution, micro-transport property, and 
8 9 

relaxation ' . These groups are in their increasing order 

of incoherence. The relaxation process, as the most ran- 
dom group, provides a memory-loss by path perturbations. 

Thus the closure lies in the memory-loss that is neces- 
sary for the transport property to reach its equilibrium, 
rather than- in closing the Reynolds stresses or other 
higher order transport functions . 

Our first task is to derive a macro-kinetic equation, 
not through the closure of the Bogoliubov hierarchy, but 
by using a scaling procedure. The result is a macro-kinetic 
equation with a diffusivity serving as an integral operator, 
so that the memory can be extended to the distribution 
function that follows. Consequently, the kinetic equation 
takes the form of a mixture of the Fokker-Planck differ- 
ential equation and the Boltzmann integral equation (Sec- 
tion II) . 
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The kinetic equation of one-point macro-distribution 
can be transformed into a continuum representation by 
means of the moment method, and subsequently, upon multi- 
plication by a macro-velocity, can derive an energy equa- 
tion for the spectral balance of turbulence, bypassing 
the need for the two-point distribution function. Since 
the diffusivity acts as an operator, the distribution 
function that follows takes a Lagrangian form, and, to- 
gether with the macro-velocity, will form a Lagrangian 
correlation. In this way, a transfer function is obtain- 
ed, describing the mode-coupling across the spectrum, 
and is founded on two Lagrangian correlations: one from 

c 

the field fluctuations for forming the diffusivity, and 
the other from the macro-distribution for describing the 
"particle-fluid" interaction. Since each Lagrangian cor- 
relation presents a memory loss, the transfer function 
possesses two memory-loss functions, the analysis of 
which must be based on the Lagrangian-Eulerian transform- 
ation. The said transformation is in the phase-space and 
differs from an analogous transformation 10-12 for dif- 
fusion in the configuration space (Section III). The com- 
petitive interplay of the two memory-loss functions pro- 
vides us with a means of closure for the derivation of 
the transfer function in the form of a gradient transport 
or in the form of a collisionless damping without a macro- 
scopic gradient (Section IV). In this manner, we are able 
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to isolate an eddy viscosity and a rate of collisionless 
damping (Section V) . 

Finally, we test the applicability of our kinetic 

method to the inertia turbulence. We derive analytically 

13 14 

the Kolmogoroff law ' and its numerical coefficient 
(Section VI) . 

II. THE SCALED KINETIC EQUATION OF TURBULENCE 
The hydrodynamical equations 

(V + u- V - vV 2 )u = E, V • u = 0 , (1) 

which describe the micro-dynamical state of turbulence, 
govern the total velocity u as driven by a field E. This 
field may consist of the pressure gradient, the buoyancy, 
or any other fluctuating forces, including the random 
noise. 

By the Reynolds decomposition, we can separate the 
total field 

E = E + E (2a) 

^ v v 

into a mean field <E> s E and a fluctuation E. In an analo- 
gous way, we can decompose the field fluctuation 

E = E (o) + E' (2b) 
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into a macro-field E (o * and a micro-field E' / with un- 
equal correlation times 

t { o) > t' , (3) 

c c 

in the increasing order of incoherence. The scaling can 
be made by means of the operators 

A = <•••>, A (o) , A' , (4a) 

which select the deterministic "global average", the macro- 
and the micro-fluctuations, respectively. The macro-fluc- 
tuation evolves in a medium possessing a transport prop- 
erty that is carried by the micro-field fluctuations. 
Finally, a relaxation is needed for the transport property 
to approach its equilibrium. This is made possible by a 
memory-loss. Thus the three transport processes of evolu- 
tion, transport property, and relaxation form our main 
framework for the closure of turbulence. For the sake of 
convenience, we introduce an operator 

A = A + A (o) , (4b) 

O 

to represent an "accumalted macro-group" of fluctuating 
field 

E = E + E (o) , (4c) 

O A* V 


and path 
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(o) 


(4d) 


z - i i + a 

o 

/V 

It is to be noted that the averages by A, A q are en- 
semble averages as specified by appropriate distribution 
functions. The groups as scaled by A, A*°^A' , only dis- 
tinguish themselves by their degree of incoherence, and 
therefore may overlap in their spectral distributions 
with respect to frequency and wavenumber. 

The deterministic functions, which result from the 
global averages of the scaled fluctuating functions, inhe- 
rit the same "rank" notations 

( )° , ( )' , ( ) 

o 

of the originating scaled fluctuating functions. 

The hydrodynamical system of equations describing the 
micro-dynamical state of turbulence is usually complicated 
by the velocity function and the driving force which make 
the equations nonlinear and inhomogeneous . In order to 
avoid these difficulties, we adopt the kinetic approach. 

This transforms the velocity function into an independent 
variable, and incorporates the driving force into the ad- 
vection in the phase space. This procedure renders the 
equations homogeneous and eliminates most nonlinearities, 
except the one representing the coupling between the driving 
field and the distribution function in the so-called "wave- 
particle" interaction. This last nonlinearity does not need 
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I 



ii n 


our immediate attention. We may postpone after the closure 
and solve by an equation of state relating the fluctua- 
tions of field and velocity. These advantages will pro- 
vide us with a simpler statistical basis of closure. 


As the point of departure, we write the following 


8 

kinetic equation : 




( 5 ) 


— o (o) 

where f(t,x,v) and + (t,x,v) are the mean and the macro- 

V * I v 

distributions, respectively, and 


l/l >t , Y* t/tx / "b s "b/W • 


The eddy diffusivity is 

/ E(t*J U(t,t -t) Efr-vjY 

* -V V v ' 

or, in the simplified form: 



V-pr Ejt-r)) , 


(6a) 


(6b) 


where the Lagrangian representation of the micro-field is 
simply written as 


= E'[t-T.x(t-V] ■ 
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and the Eulerian field at time t is written as 

= B'(t ) 

s\j ^ 

since x(t) = x in E / (t / x(t)) = E ; (t) . The propagator, or the 

'V -V 

Lagrangian evolution operator U(t,t-T) describes the evolu- 
tion of the perturbed trajectory. The diffusivity may serve 
the role of an integral operator if it is written as D'{}. 

III. THE LAGRANGIAN- EULERIAN TRANSFORMATION IN THE PHASE SPACE. 

THE LOSS AND CUTOFF OF THE MEMORY 

By means of the Fourier transformation in k-space, 

-v 

the Lagrangian correlation can be written in the following 
form: 



It contains the spectral function 


* y ( e.'(ic) B'(- k’)y 

X ~ \ <v V 

giving the mean intensity 

fjLk’<^'(kV) = frw> - 

& ** ** 

where 

y * tyxf 
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is a factor of Fourier truncation within a length interval 
2X in d dimensions. The Fourier transformation has the 
Fourier kernels 


4i (-T,k"j s -t kT , b (-z, ky) s -e 

CO /v ' ^ v 


- Ck" v' T 


( 8 ) 


The frequency w k = u> k (k") is a function of k" as determined 
by a certain dispersion relation, for example, reflecting 
the interaction between internal gravity wave and turbu- 
lence 15-17 . The perturbed path l( -t) forms a fluctuating 

j. k ;/ • ( — t ) 

orbit function e ~ -v . Since we can decompose the path 
l into groups 

, 


we find the corresponding "scaled orbit functions", or 
briefly, the "orbit components": 




(9) 


On a separate occasion, we have developed a theory of 

transition probability for the orbital motion and deter- 

18—20 

mined the following correlation function to be : 







( 10 ) 


The first two factors on the right-hand side represent the 
factorization when the fluctuating functions 
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and 


( 11 ) 


/ 


4 



are mutually independent, and h^ is a correction factor 
when the hypothesis of independence does not hold. 


With the substitution of (10) into (6), we find that 
the Lagrangian correlation 


w'(c) x)e (t-t-jy = fdk" C}’(kt Mfrky) ( 12 ) 

■v. ^ -v ' J ~ V -v ~ ~ "V 

has, in the Fourier integral, a product of the spectral 
distribution <s / ( k")> with the orbit function 

« -V 

M = L'L & 4^ 0) t' . (13a) 

to v £ 


The components may form the groups 


and < 13b > 

representing the streaming by v,u,E and the memory-loss 

v "V -V 

by the macro- and micro-fluctuations, respectively. For 
the sake of convenience, we write the micro -component 

< 13 =1 

that incorporates the inter-dependence between the fluc- 
tuating intensity s' and the path fluctuations S/; also 

* -V 

we write the accumulated macro-component as 


B-ll 



(13d) 


\ s ^ £ \ 

The function M will be called the "memory-loss function". 

In particular, if the streaming predominates, we can 
degenerate (13a) into 

M = , with M = I , 

so that the Lagrangian correlation (12) is reduced to the 
Eulerian correlation by the relation: 

W L (x) = W E (zJ = (e '(t,x) *-%)} 

u A/ v 

2 * 

J - w ~ (14) 

This is the well-known relation under the hypothesis of 
"frozen turbulence", and is valid in weak turbulence. In 
the general case where the turbulence is not frozen, we 
have the relationship 

Vjfy = W/jM(v)j ns) 

between the two correlations by using W'{} as an integral 

E 

operator. Note that 

W L '(v) 4 

% z 

since M < 1 . 
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In the derivation of a transport property or a trans- 
port function in equilibrium upon the x -integration of a 
Lagrangian correlation, it is obvious that the memory 
should not perpetuate indefinitely but should ultimately 
be cut off, as a requirement of the closure. We see that 
the memory function M consists of two components: The 
macro-component h (o) incorporates the strength and the de- 
cay of the memory. The micro-component h^ , as belonging to 
the most random group at the tail of the memory-chain, has 
the function of cutting off the memory as the result of 
the diffusion by the micro-fluctuations. 


Recall that the diffusions by u- and E-fluctuations 
are governed by the dynamical equations of the orbital 
motions 


and 


A 



4t 


dX f ol u *2 

— n_ =■ M. , - = E 

*t ~ 4t 


(16a) 


(16b) 


respectively . 

We perform the scalings and the ensemble averages by 

18 19 

means of the scaled probabilities of transition ' 

The dynamical system (16b) gives the paths 


€'C-r) = - 




Z * 


Z I 


E (o V 


(17a) 
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and a macro-variance 


belonging to the unmatured diffusion, or small-time diffu- 
sion. 

The mi:ro-variance can be written in the form 


(%(-T) = 2-K't , (18a) 

V /V A? 

belonging to the matured diffusion, or large-time diffu- 
sion with a transport coefficient 



(18b) 


Two forms of propagators U(t,t-t # ) may be visualized, de- 
pending on the orbit evolution, whether it be perturbed by 
the u-f luctuations or by the E-fluctuations, giving an ed- 

A# A, 

dy diffusivity or an eddy viscosity, respectively. Alter- 
natively, the micro -variance in the diffusion by E-fluctua- 
tions following (16b) is expected to deliver an eddy vis- 
cosity too, in an indirect way through D' and after the 

aJ 

closure. For this reason, we shall assume a memory cut-off 

and a closure by (18a). A rough evaluation of K' could be 

at 

made by relating it to the noise spectrum of the Brownian 
problem, or to the diffusion model (16a). But from a real- 
istic viewpoint, K' must be related, in a self-consistent 



way, to that same eddy viscosity which governs the mode 
transfer by small scale fluctuations (Section V) . We adopt 
the latter viewpoint. 


The variances (17b) and (18a) determine the orbit com- 


ponents : 

i (?) -d C ° J ,f -cl' 

^ = -L f (19a) 


with 


4' - -irk* A' (-*!)■ 


(19b) 


It is to be mentioned that, unlike h' which depends on 
the variance <£/ ( - x) £' ( - x)> of one path l 1 ( - x) , the 

" V M 

component 

s t' (20) 

contains a function E,' that depends on the micro-fluctua- 
tions of two paths, one is V ( - x) and the other is 
i' ( - x ) carrying the field intensity s' , so that 

S A/ 

/V 

? <E '(o)E'(-c})> . ,21, 

* V <v 

For small x, the two paths have still to grow, so that 
is small. For large x, the two paths are too far apart and 


B-15 



lose their correlation, so that tends again to zero. 
It is for the intermediate x that E,' has a finite value, 
so that (20) may even become negative, yielding a nega- 
tive portion of the Lagrangian correlation (15). 


Now we integrate the Lagrangian correlation with res- 
pect to x, to find the diffusivity in the phase space: 



j j-c Uk" ■ 

V) yJ v ^ v 


( 22 ) 


This expression can also be used as an integral operator 
for integrations with respect to x and k w . The memory- 
loss function M, as forming a part of M, governs the loss 
of memory for the diffusivity to approach its equilibrium. 
The time integral 


[ clx M(-Xk\) 


will be called the relaxation time. 


(23) 


IV. CASCADE TRANSFER WITH DOUBLE MEMORY-LOSS 

We take the moment of the kinetic equation (5) to 
obtain the following moment equation: 


^ (°) _/_ (o) (oj - (©; Co; (©; . 2. {o) -T-i 

\ *V + y(u~ U. -M u- Cc- J-W U. / =■ T 

* 6 v * </ c j- * y J i ■, 


-ty 


(24) 


with 
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( 25 ) 


M- 

Upon multiplying (24) by u.f 0 ^ and averaging, we obtain 
the energy equation: 

= P ( °V c w - T‘° ; - p-l”. < 

where the transport functions are as follows: 


production 

coupling 

transfer 

dissipation 

flux in energy 
transport 



(27) 


The momentum equation and the energy equation are not 
purely hydrodynamical , since they contain the distribution 
function f*°* through and T*°*. The function rep- 

resents a collisionless interaction, and the transfer func- 
tion T*°* represents a transfer between the two portions of 
the spectrum of macro- and micro-energies, respectively. 
This transfer originates from the "particle-fluid” inter- 
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action, as shown by the new Lagrangian correlation 


^ 0 -(-r,yJ= (28) 

between the fluid velocity u^°* and the distribution f*°* 
of the "particles". Thus the transfer function is more 
specifically 


J*>> 

V V 

< 29 ’ 

from (25), (27) and (28). 

By applying the rule (12), we can transform the new 
Lagrangian correlation into 


with 



(30a) 



(30b) 


and get the transfer function in the form: 


T 1 = -Uk'Uv v. ■& • yl a (31 

d ” J ~ v <v ' 


This form can be made more explicit by inserting the ex- 
pression (22) for D* , as follows *• 
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(32) 


r‘- Nji-xk-t •) 


Here we have introduced for the sake of convenience: 


- < 33ai 

or, more briefly. 


V J ~ * /< ~ ' * * u ~ J f 


(33b) 


where the dependence of M(k")/ M(k ; ) on -t,v, and of 

v ' «. 

, <p!°* (k') on v may be considered to be understood. 

By expanding the multiple derivatives 3^,3 r , we have: 


2 fo i [«&*; fM ' 1 $y+M(ykky jlffip 

+ Shvl ■ 


(34) 


with 


V V ^ A 


+ \Mtki Mkn$k'ji-M(ky \kki)f-w) 

y ~ ^ ^ * ~ 4 * 


h(iy 

- 4 ~ 


(35) 


The second term on the right-hand side of (34) gives, 
after two successive integrations by parts: 
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( 36 ) 


Jh /% if )\ 

with the obvious condition 

= ° for I V’ I ■» OO , 

A/ 


so that, after a substitution, (34) becomes 



As far as their real values are concerned, 

N 'rt&/kV *1%'), M(k") and M(k') are even functions of 

k" and k # , so that the derivatives 

7> M(k°J and ~b M(k') 

will generate odd functions, and will not contribute to 
the integrations with respect to k" and k' in the infinite 
domain. Therefore they can be deleted from the expression 
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for the purpose of calculating T 


from (32) . 


of N 


(o) 


3 r 

The result 


(o) 



( 39 ) 


is in the form of the moment of tp|°* as driven by the 

differentials 3.3 of the orbit functions. Hence the trans- 

3 r 

fer function (32) takes the form 


« f(k'i pry . W) 

and is seen to be controlled by two memory-loss functions 
M(k ,; ) and M(k') as imbedded in M(k") and M(k') . 


V. EDDY VISCOSITY AND COLLISIONLESS DAMPING 

It is to be noted that the v-dependent orbit functions 
have their slowly varying component h* and their rapidly 
varying component h v> Thus the differential 3M will be 
mainly contributed by the differential 3h^ of the rapidly 
varying component in the manner 


~2>M = -Lkr M 


( 41 ) 


We notice that the two orbit functions 


M(~T,k'v) and A/ (~h k' *) 
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compete in their roles of loss of memory for the approach 
of the transfer function (40) to equilibrium, under the 
scaling conditions 


0 4 i' £ k £ fc'V Co 


(42) 


Thus it suffices to select the function that is most ef- 
fective in this role. On this basis, we shall calculate 
the asymptotic expressions by considering small and large 
scales. Since the path amplitudes and the trans- 
port coefficients ,K' are separated by k, the small 

3 * 3 * 

and large values of k correspond to : 


P > V 


(43a) 


and 


r< X' , 


(43b) 


respectively. 


Case (a) . Large k 


For large k, or for even larger k" , such that 


k" » k 


(44) 


in conjunction with (42) and (43a) , we can attribute the 
dominant role of memory-loss to the function 
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(45) 


fo(-x,k ") % s ') 


which is independent of v, while leaving the role of 
streaming to the v-dependent orbit-function 

A Itetl * ■iJ-r.k'J l v (-T.k',v) . < 46 > 

The independence of (45) on v is seen from (18) and (19), 
after having assumed = 1 or 5'sO . This approximation 
is legitimate, because in the integration (40) with re- 
spect to x, the role of hi. in M(-T,k",v) is seen from 

c, «■ » 

(20) and (21) to be taken over by h 1 and h Q for small and 
large t, respectively. Hence we reduce (39) and (40) to 


N 


(o/ 






(47) 


and 


T 


*) 


k. 


•/«* (^i> Jiy K- rfcMW;), ^ 


respectively. Here 


) i(-T. t; v) U-- 

*o 

is a modulation function having a streaming by h v and a 
memory-loss by h Q h' • 


y)tjrX.k-)L'(^X) 

** ° -v V 


(49) 


Cf(k,k 
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In. strong turbulence, the memory-loss predominates 
over the streaming, and renders the modulation function 
to become independent of v as in 

A* 

Q( k ') -fd t t 2 Ij-t-x) • < 50 > 

- u 0 * ° « 

Thus we obtain an eddy viscosity 


and transform (4s) into 



K.' fa' M fa* 

rJ * i /•*/■*•* 


9 %» 

* <V -V 


The integrations herein involved are: 


(51) 


(52) 


A ^ -A * * 

IV v 

- a <w <n-;> 

~ -v 

jfa A% : K'<W>' 


(53a) 


(53b) 


and 


A A * 


-<: 


y * A t 


(54) 
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Hence, by substitution into (52), we obtain the transfer 
function 


Y° J K' 


(55) 


in the form of the product of the eddy viscosity K' r by 
the vorticity function . Here, by isolating an eddy 

viscosity, we have found a cascade transfer of the gradi- 
ent type . 


Case (b) . Small k 

The small values of k imply an even smaller k # . This 
extremely small k' , as compared to k" , suppresses all 
roles of M(-T,k',v) before M(-x,kf',v)r rendering 

^ -V ** AS 



(56) 


Thus we transform (39) into 


N* J 


tnidJUL k 


0 V ^ « V V " 


(57) 


The factor k"k" emphasizes large values of k" so that 
we can write 

K *) = A^ k v \ (58) 

as an approximation. When this function is substituted 
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into (57) and subsequently integrated with respect to t. 


following (40) , the contribution from h v as a streaming 
is again negligible as compared to the memory-loss for 
large k #/ in strong turbulence. It follows a modulation 
function 

Q(k ") = f Jz T 2 ft) (-% k *) 

~ J o 

= ft* rUj-T,kU(-T,k'jir-T,k'') 

° (59 

for the control of the transfer function (40), yielding: 



Here the two integrals with respect to k' and k" become 
separated into an energy 

(& l \ - [&[***< ¥%;*) (61a) 

by (56) and (30b) , and a frequency 

nr . /*• k" Ujrh $(r) . mb, 

v q / v -v 

called the rate of damping in the collisionless damping 
of the energy. Hence the transfer function becomes: 



( 62 ) 
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Since the transfer of energy is not of the gradient 
type, we cannot isolate an eddy vicosity as was the case 
of large k in (55) . 

We may consider 

} and 

as integral operators which emphasize large and small 
values of k by their respective integrations 



and write 


k :\ } 

V*' 1 



(64) 


In this manner, the two transfer functions (55) and (62) 


become 


-f) 




= K' 


r 




and 


T < 0) 



(65) 


( 66 ) 


For large and small k, the operators (63) pick up 



respectively, from the sum k£k'_. + k/'k " as appearing in 

V v ^ -w 

of the general transfer function (40) . Thus the formula 
(40) is a general transfer function with the special forms 
(65) and (66). 

VI. DERIVATION OF THE K0LM0G0R0FF LAW FROM THE SCALED KINETIC 
EQUATION 

Although the Kolmogoroff law 13 of the inertia turbu- 
lence has been derived by many investigators using diverse 

2 1 

methods as seen in a recent review of the diagram tech- 
22 

nique in conjunction with the direct interaction approx- 
2 3 

imation , it is for the first time that the analytical 
derivation of the Kolmogoroff law and its numerical coef- 
ficient is obtained from the scaled kinetic equation of 
turbulence. Our kinetic method is generally valid for in- 
ertia and non-inertia turbulence. 

We assume an isotropic and homogeneous turbulence, 
and that the field 

f ,67) 

comes from the pressure gradient Vp in an incompressible 
fluid of constant density p. The shear and the buoyancy 
are absent. The condition of incompressibility, as applied 
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to the Navier-Stokes equation (1) gives the relation 

VE =■ V7 - ua < 68) 

V v V y % 

between E and u, and therefore also between their spectral 
distributions S(k) and F(k). We have 

t <C>=j!kS(k). t<^>-f7kr(k), 

o o 

and find the relations 


f pk n(k) , 

(69) 


$(*,;= < F(k) T?>(k ) 

n Lk) = f &- z s(k). 


where R (o ^® R^° ^ is the vorticity function, from (54), and 
the numerical coefficient is determined td be 19,20 



The details of the calculation are omitted. 

In the inertia subrange, the energy balance in the 
spectral form is 13 ' 14 

= e - (72) 

describing a constant cascade of energy transfer across the 
spectrum at the rate of energy dissipation e. Now the trans- 
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fer function, as found in (55) , can be written in the iso- 
tropic form: 




(73) 


where 

k'- fj^k‘ $(k"J Cf ( k"J (7 

is the trace of the eddy viscosity tensor (51), and the 
modulation function is given by (50), with ( 0 ^ = 0. 

From (9), (18a) and (19), we rewrite the modulation 
function (50) in the form: 


^ f(k, K ;-/3t T Z J 


(75) 


It contains the decay of the memory at a rate 

(?) 


^ = r_L k " 2 (e^S 
_ 2< / - _ 


(76a) 


and the cut-off of the memory at a time to' 1 such that 


< = k° L K\ 


(76b) 


The integral is evaluated by an interpolation, giving a 
modulation function 




- ft 


'«k } 3 


(77) 
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Hence we determine the eddy viscosity by the following in- 
tegral equation: 

K' • f fd,k" S(k,“) Q [*" . < 78 » 

When the spectrum S(k) is converted into the spectrum 
F(k) by the relation (70), we find the solution of (73), 
giving the spectra 

F(fe ) - A k ^ with A = /• 7 
S(k) = j A z L* 1 ' 3 k~ 1/3 

- n(k) - f A z f x k- 7 *. 

(79) 

of velocity, field and pressure fluctuations, respectively. 
The details of calculations are omitted. 

VII. DISCUSSION 

The micro-dynamical state of turbulence can be de- 
scribed by a nonlinear and stochastic hydrodynamical system 
which may be homogeneous or inhomogeneous. The Navier-Stokes 
system of equations is inhomogeneous, as in (1), or briefly 
as 
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(^+L)ul = E, Z^ = °, With L <801 


The inhomogeneous system in the physical space may be con- 
verted into a homogeneous equation in the phase space as 


N(t,*y)=0, 

V V 


with v/y-Vt 7 ■+■ £ .<) 


(81) 


called the "master equation". It is equivalent to the in- 
homogeneous system (80) , provided 


N(t,X,v) rn j> S [v - ufcxj] . 

V -w V V 

Other homogeneous systems may be: 

(i) the Burgers equation 

+ L) u = o, with L . u _ v 2l , 

(ii) the Korteweg-DeVries equation 

: 

(W) u. = with L - LlA. +/ xAl_ 

N r 7 -^>X 

(iii) the equation of geostrophic turbulence 

(\ i- L) V<U = 0 , with Lh a- \7- V F 2 . 


(82) 


• (83a) 


(83b) 


(83c) 
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From the conceptual viewpoint, being given the micro- 
dynamical state of turbulence, the question may arise as 
to what treatment one should choose: the kinetic treatment 
or the hydrodynamic treatment. Since a homogeneous partial 
differential equation can describe the exact orbit dynamics 
by the differential perturbation operator L alone, there 
is an essential advantage to keep a homogeneous system, 
either in the hydrodynamic representation (83) or in the 
kinetic representation (81). For this reason, we have 
chosen the kinetic method which, moreover, supresses all 
nonlinearities caused by the velocity since the velocity 
is an independent variable here. The only nonlinearity 
that subsists is from E- 3N in (81). This entails a Rey- 
nolds stress in the phase space, the evolution of which 
requires the knowledge of the Lagrangian correlation of 
E-field fluctuations and a theory of transport. The ap- 
proach to equilibrium and the loss of memory constitute 
the basis for the closure of turbulence. 

There are two Lagrangian correlations involved: one 
is the correlation of E-field fluctuations as stated above, 
and the other is the correlation between the velocity fluc- 
tuations and the fluctuating distribution of velocities. 

The latter arises because our kinetic equation (5) of tur- 
bulence has a long memory. Since each Lagrangian correla- 
tion entails a memory, our kinetic theory of transport 
finds two memory-loss functions. Their competitive roles 
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I 


and their cut-off are the main issues of our analysis. 

The difficulties are minimized by introducing a scaling 
procedure which dispenses from our involvement with the 
two-point distribution function. Then the scaled proba- 
bilities of transition determine the scaled orbit func- 
tions. The scaling permits the distinction between the 
macro- and micro-variances of path fluctuations which are 
associated with the unmatured (i.e. small time) and the 
matured (i.e. large time) properties of transport. Other- 
wise, the unsealed variance would appear as an external 
parameter. 

Our results indicate that the transfer function takes 
a gradient form for small scale turbulence, isolating an 
eddy viscosity, and a form of collisionless damping with- 
out gradient for large scale turbulence, isolating a rate 
of collisionless damping. 

We can write the equation of energy (26) in the form: 

\j k ^F(k') - p l ° ; - tr- t v ~s j . <84, 

with the transport functions: production , loss B*°^ , 

transfer T*°^ and dissipation e^°^. In statistical equili- 
brium the left-hand side is independent of k, so that, 
after a differentiation with respect to k, we have the 
spectral balance 
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(85) 



-(P - 3 ') 


for scales larger than that of the inertia turbulence. 
Note that e ceases to be the governing parameter here. 
This means that 



> o 

in a net production 

(86a) 

f(°j 

< o 

in a net loss. 

(86b) 


The change of signs in (86a) and (86b) refers to a direct 
cascade (i.e. a cascade transfer toward large wavenumbers) 
and to an inverse cascade (i.e. a cascade transfer toward 
small wavenumbers) , respectively; The transfer function in 
the form of collisionless damping is preferred for the de- 
scription of the inverse cascade. The arbitrary use of the 
transfer function in the form of gradient transport would 
result in a negative spectrum. 
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CHAPTER C: QUARTER THREE 


Group-Scaling Theory for the Enstrophy Transfer 
in Two-Dimensional Turbulence 


C.M. Tchen 

The Graduate Center and The City College of 
The City University of New York, N.Y. 10031 


ABSTRACT 

Since the detailed interactions between the indi- 
vidual modes contain too many minute details unneces- 
sary for a statistical theory, we scale the modes into 
groups which can more easily decipher the governing 
transport processes and the statistical characteristics 
An equation of vorticity transport describes the micro- 
dynamical state of two-dimensional, isotropic and homo- 
geneous, geostrophic turbulence. By the group-scaling, 
we derive the equation of evolution of the macro-vorti- 
city in the form of the Fokker-Planck equation with 
memory. The eddy diffusivity and the enstrophy trans- 
fer relax to equilibrium through the memory-loss by tur 
bulent dispersions. The memory-loss is analyzed by de- 
veloping a theory of probability of retrograde transi- 
tion in the scaled form. We find that the eddy dif- 
fusivity contains two time scales, characteristical of 
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the small-time scattering and the large-time diffusion. 

In the special cases, our formula can be degenerated to 
an empirical form suggested by Heisenberg when the scat- 
tering prevails, and to a form derived earlier by Tchen 
when the diffusion is dominant. In the inertia subrange, 

as governed by the enstrophy cascade, the spectral law 
2/3 - 3 

F = Ce _ k is derived with a numerical coefficient 

u C 

C = 2.59. 
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I. INTRODUCTION 


The micro-dynamical state of turbulence can be de- 
scribed by the equations in the following forms: 

(a) a non-conservative (or inhomogeneous) equation 
in the physical space, 

(b) a conservative (or homogeneous) equation in the 
physical space ,, 

(c) a conservative equation in the phase space. 

The form (a) contains a stirring force, which may be 
the pressure gradient, the buoyancy, or the random noise. 
The Navier-Stokes equation falls under the form (a) . The 
two-dimensional geostrophic turbulence also belong to 
the form (a) , if the turbulence is driven by a random 
vorticity source, and will fall into the form (b) if such 
a source is absent. The one-dimensional Burgers equation 
and the Korteweg-DeVries equation also falls into the 
form (b) . 

It is to be noted that the non-conservative equation 
(a) in the physical space can be converted into a con- 
servative equation (c) in the phase space by regarding 
the stirring field as an advection. Therefore the prob- 
lems of the conservative system should be the main ob- 
jective of a statistical theory. With this emphasis in 
mind, we develop a statistical theory, using the scaling 
procedure to study the mode-couplings and the transport 
property, and relying upon the memory-loss to obtain the 
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closure of turbulence. A kinetic theory 1 on the basis 
(c) has been developed by Tchen for the Navier-Stokes 
turbulence (a) . Presently, we apply the scaling method 
to the two-dimensional, isotropic and homogeneous, geo- 
strophic turbulence. 

In view of the many applications to the meso-scale 

flows in the atmosphere and the availability of numeri- 

2-7 

cal computations , the problem of geostrophic turbu- 

8—10 

lence has attracted many investigators , see a review 
by Rhines 11, The geostrophic turbulence is characterized 
by an enstrophy transfer from large scales toward small 
scales, yielding a spectrum of velocity fluctuations 

_ 3 

following the k law. Most theories analyzed the triad 
12—18 

interactions 

The present theory simplifies the interactions by 
scaling a fluctuation into three groups, representing the 
evolution, the eddy transport and the relaxation. Our 
purpose is to derive analytically the enstrophy transfer 
function, the eddy transport coefficient and the spectral 
distribution. Upon a generalization to a kinetic method 
(c), as is equivalent to the non-conservative micro- 
dynamics (a), we expect to find: 

(i) the spectral law k _5,/3 in a direct cascade of 
energy transfer at large k, 

-5/3 

(ii) the spectral law k in, an inverse cascade 

of energy transfer at small k. 
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(iii) the spectral law k from an enstrophy trans- 
fer in a geostrophic turbulence that is driven by ran- 
domly distributed vorticity sources. These results will 
be reported on a separate occasion. 

II. SCALED VORTICITY EQUATION OF TURBULENCE 

The two-dimensional incompressible flow in the hori- 
zontal plane is governed by the following equations of 
motion : 


(\ + l ) ? (%*) = o ,1a) 

and 

v- u - o ( 1 b) 

-v" «■» 

Here u is the velocity, 

£ = Vx Lt - ° 'K ) ( 2 ) 

is the scalar vorticity, 

L = U V -V/ (3) 

v -v- 

is the differential operator, and v is the kinematic vis- 
cosity. By the condition (lb), the pressure gradient, 
the vertical buoyancy force, and the vertical Coriolis 
parameter are eliminated. The system (1) describes the 
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micro-dynamical state of geostrophic turbulence. 

The scaling may be made by means of the operator 


A~<> 


(4a) 


of global average and the fluctuation operator 


A * I- A 


(4b) 


giving the deviation from the global average. The fluc- 
tuation operator 


A = A (o) t A 


( 4c ) 


is subsequently divided into the operator of macro- 

fluctuation and the operator A' of micro-fluctuation. 

For the sake of convenience, we introduce an operator 


/) = A+A ( °>= I- A' 


<4d) 


giving an accumulated macro-group. Thus we have: 


A It = U , 


U + u- 


(5a) 


and 




(5b) 
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The macro-fluctuation and the micro-fluctuation can in- 


stantaneously overlap in their wavenumber components, but 
their statistical functions are distinctly separated by a 
wavenumber variable k, and by their durations of correla- 
tion 



>^c 


( 6 ) 


in the increasing degrees of incoherence. 

By applying A (o) and A' to (la), we obtain the follow- 
ing scaled system: 

(S t * A ( ° J l) = - V- /V? ' (7a) 

(■^ + L ) K' = - u'F V- • (7b) 

By a formal integration of (7b) and a substitution 
into (7a) , we eliminate £' to obtain the equation of 
evolution of the macro-vorticity : 

v-K - i« 

where 

K' =j t *?T (lift x) Vfct-T) (9) 

is the eddy diffusivity, and U(t,t-x) is the propagator. 
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or the operator of evolution along an orbit that is per- 
turbed according to (3) . Without altering the value of 
the integral, the upper limit may be put to infinity, by 
the scale separation (6), giving a matured (i.e. large- 
time) diffusivity. It is to be noted that the diffusi- 
vity enters as an integral operator k ' () in (8), so that 
this scaled vorticity equation takes the form of an integro- 
differential equation, and not the customary Fokker-Planck 
differential equation. 

When we multiply the vorticity equation (8) by 
and average, we get 

t V? n ) * £?((£ +<i o> )X. (0)2 > 

Here we have retained the kinematic viscosity in the dif- 
ferential operator L. 

In a homogeneous turbulence, the second term on the 
left-hand side vanishes, and the last term on the right- 
hand side becomes 

t($ (o) =-e ( ° ; ^ith . (11a) 
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In an analogous way, we can write 


= -K<r >l > ■ (i ,b > 

and find the first term on the right-hand side of (10) to 
be 

- <s v-K'i Z s’"}) =K'{' (ftti ■ if{t- z P}* T*’- die) 

Here the diffusivity remains to be an integral operator 
with its time integration extended to the Lagrangian cor- 
relation 

(fttl ■ (lid) 

^ V ^ 

The diffusivity is written in the form of a trace. The 
function is the rate of molecular dissipation of the 

enstrophy i<5^°* 2 >. The function T^°* represents the mode 
coupling and is called the "enstrophy transfer". 

In terms of the dissipation function and the transfer 
function, the equation of spectral balance (10) is reduced 
to the form: 

It represents the time decay of the mean vorticity or 
enstrophy 
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(13) 


±(^ >Z ) = J k dk' F ? (k1. 


which has a spectral distribution (k) of £-f luctuations . 


III. TRANSPORT PHENOMENA IN TURBULENCE 

A. Lagrangian-Eulerian transformation in the scaled form 

The Lagrangian velocity follows a trajectory that is 
governed by the dynamical equation: 


)_ 


- «.(t,) 


(14a) 


or the orbit equation 

jM-*) 




= - , 


(14b) 


where the perturbed path 


\(-t) = J(-z )+JU-t) 


(14o) 


consists of a mean path Jt(-x) and a path fluctuation 
£(-t). In the following calculations of the path pertur- 
bations, we shall neglect the effect of the kinematic 
viscosity v as compared to the larger effect from the 
eddy motions . 

By means of the Fourier transformation, we can write 
the Lagrangian and Eulerian correlations in the following 
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forms : 


K'(i) * fTt *; vfcte) 

~ - 

<Cc-t,*'y (' o ' m *' * '' fvJ ) n51 

^ <w ^ 

•m 

and 

V.6v> s' * ) l*(t-T, X-Jj) 

^ v - — 

= Uk"ij-z,k"j (^jkv'yU-i.ki, !(-T,k ¥ ht k ~ ^d6) 

* ~ - 

respectively. Here sMk") is a fluctuating intensity in 

2r 

the k-space, such that 

(<#•/> = / (u-(f/ ’iC-kiy ,i7a > 

is the spectral function, and 

4. ^M<V(r/> = (<{(ubu(u)\.zpk- Fjki <17b) 

* " % 

is the mean intensity in the x-space. The coefficient 

is the factor of truncation, when the Fourier transforma- 
tion is truncated within a length interval 2X in two 
dimensions, i.e. d = 2. The Fourier transformation has a 
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Fourier kernel 


• n8) 

with a frequency w = w (k") that arises from a dispersion 

JC Jv 

1 9 

relation, e.g. as relating to an internal gravity wave. 

The integrations in (15)- (17) and in the following 
are understood to extend from - 00 to 00 . 

The correlation in the integrand of (15) can be 
written as 



consisting of a mutually independent part that can be 
factorized as 



"V V / 




£ 



A 




(20a) 


and a correction 



(20b) 


C-12 



accounting for the mutual dependence. We find that this 
correction is not too important for the small-time scat- 
tering and the large-time diffusion that are of interest 
in our problem. The mutually independent part contains 
an average 

A A, -e * * (2i) 

0 


that is calculated by the probabilities of transition in 
the scalings by A q and A. The derivation of the proba- 
bility functions and the calculations of the average are 
presented in Appendix A, finding 

AA 0 e Ck -"- H = l(-xkU°Uk')l‘(-Tk1. < 22 

V ~ 


On the right-hand side appear the components of the orbit 
functions, or "orbit components". The mean path gives the 
component 



-£ 


ck'' i(-*) 


(23a) 


The fluctuations give the components 

'Cf-'Cr/ = (f- ^ (_ c^xj ^ 


C-13 


I 



in the unmatured and matured dispersions, respectively, 
of tilne scales m^°^ and , such that 



<o' = k"k.".- K' s- k" l K' • < 25b > 

A* M V 

V 

The details of the calculation are given in (A20)-(A2 4) 
of Appendix » The function F u (k") is the spectral dis- 
tribution of velocity fluctuations, as described by (A19) 
in Appendix . 

If we lump the orbit components caused by the fluc- 
tuations into the function 

m * h CoJ r , 


called the "memory-loss function", we can write the Lag- 
rangian-Eulerian relationship in the form: 

k'c.j = u/ £ y Mh,ki} • (27) 

This is a nonlinear integral equation for W' (x), since M 

a, L 

a* 

contains through h' , a/ and (9) . Note that the Euler- 
ian correlation is governed by the mean path £(-x), while 
the Lagrangian correlation is governed by the mean path 
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and the path fluctuation £ (- t) . 


B. Relaxation time 

Upon substituting (27) into (9), we find the diffusi- 
vity tensor and its trace in the form: 


(^u^i "J) J h k 7 

-V* W .* v 

~ o 

* /4 ' (^L &7) £ ^ 7 ' 

v a/ K ^ 

K'- &/<'_ fcU." L Gr (k-J 

Jv / ^ ^ v A** <v 

= 3 -/^" 


(28a) 


(28b) 


It is seen that the diffusivity approaches its equilibrium 
in a relaxation time 


= [ <jr M(-~c, k "J 

-f*v SL*p(-y»°* z -c z -c^~c). 


or, approximately. 


-I 


Q-Jk'U . 


(29a) 


(29b) 


The approximation is based upon the interpolation. 
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C. Diffusivity 

The equation (28b) which determines the diffusivity 
J^' is in reality an integral equation since the relaxation 
time G (k") is also a function of /<;' through to', by (25b). 

As an approximation, we calculate the asymptotic 
values of /C ; by considering the following two special 
cases : 

case (i) If G_ is controlled by n/°^ , we find 



with = VTt in two and three dimensions. 
t\ 

case (ii) If G^ is controlled by to' , we find 



(31) 


As an interpolation, we write: 

K'-jj-fik' Fjk"j( \±- ! k" 3 F u ti")\ t . 

k ^ 



(32) 


It is degenerated into 


K' = 

** 

• U K« 



L K' . u 

K) j 

t •%■» 

■* C °> 


( 33 ) 
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IV. DOUBLE MEMORY-LOSS 


The first memory-loss function M(-x,k") appears in the 
Lagrangian correlation (27) which produces the diffusivity 
(32). The second memory-loss function arises from the 
second Lagrangian correlation (lid), with the same dif- 
fusivity in its role of an integral operator as appearing 
in the transfer function (11c). 

In order to derive the structure of the transfer func- 
tion (11c) in isotropic turbulence, we first write the 
second correlation function as 

('fit)- 

in analogy with (27). Here 

is the spectral intensity of ^-fluctuations, such that the 
intensity 

= <r%) -ffriJ) 

can be expressed in terms of (k ; ) , by 

'(rfkk'j) = 2-fi-k ' i' ') 

* 7 l L °’ . < 36a > 


(34) 


(35a) 


C-17 



Use is made of (11a). Note that the relationship 

ty) = ± k* /; (kj (36bi 

between the two spectra follows directly from the defini- 
tion (2 ) . 

Now we apply the integral operator f\' { } to (34) and 
write the transfer function in the following form 

= K '[[ . ( 37a ) 


according to the definition (11c), or in the alternate 
form 



(37b) 


when K'{ } is written out fully by means of (28b) . It is 
seen that the transfer function approaches its equilibrium 
at a new relaxation time 

Q-(k',k") =T M(< k') M(-T, K V , (38) 

T ~ ~ J 0 

with two memories: 


M(-t, h'J h\ (-T, k "J . 

A. ** 
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V. TRANSFER FUNCTION 


By the use of (13) and (17), the transfer function 
(37b) can be rewritten as: 



where, by definition (38), the relaxation time G T (kJk M ) 
sums up the contributions from the two competing memory- 
loss functions (39) under the conditions 

k ' < k < k " (41 ) 

with the properties 

hi <*Ji cocos') « <42) 

so that the role of memory-loss is taken over by M(-T,k") , 

i . e . 


M(- T,k') " I • M3) 

Th u s we can reduce (38) to 

Q T (.Kk") s Q„(kV. < 44 > 

and separate the two integrations with respect to the wave- 
numbers in (40), obtaining the degenerate form 
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( 45 ) 


T^> -ffik' FJWW) z Ak' k'^ik'J . 

k 0 

or, more briefly, 

T t f- K ' ^ . 

in the notations (28b) and (36a) . 


(46) 


VI. INERTIA SPECTRUM 

In terms of the spectra F u (k) and F^(k), we can write 
the spectral balance (12) in the form: 

ifak ' 1 F;(k'J ■ (47) 

* o *0 

Since the differential vanishes, i.e. 


1 W = o 


in statistical equilibrium, the integral form should read 


(K + ») ‘J <lk'k' z F Jk') . 


(48) 


where 


£ 



(49) 


is a constant rate of dissipation of the enstrophy. A 
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fluctuating would enhance the intermittency which we 
shall not consider at present. 


The inertia subrange is described by the spectral 
balance 



(50) 


by omitting the effect of molecular dissipation in the 
left member of (48) . 

When we use the general formula (32) for fX , and ob- 
serve the relationship (36b) between the two fluctuations 
u and £ , we determine the spectra from (50). The results 
are: 


Fjk)~Ci^k 3 , C=z--f 

c<V\ 


(51 ) 
(52) 


VII. DISCUSSIONS AND CONCLUSIONS 

The problem of turbulence starts with the nonlinear 
equation describing the micro-dynamical state of turbu- 
lence. This equation contains too many minute details 
which are too complicated and also unnecessary for the 
statistical theory of turbulence. Past theories either 
followed the kinetic method and closed at a high order 
distribution function, or the hydrodynamical method and 
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closed at a high order correlation. They analyzed the 
detailed interactions between individual modes. The 
present theory begins with the coarse-graining procedure 
by distinguishing between the fluctuations of the macro- 
group and the fluctuations of the micro-group. It trans- 
forms the equation describing the micro-dynamical state 
into a scaled equation of turbulence, in which the macro- 
fluctuations evolve in a medium of more random eddies be- 
longing to the micro-group. In this way it analyzes the 
statistical couplings between groups of modes instead of 
the detailed couplings between individual modes. For the 
approach of the transport property to equilibrium, as a 
closure of turbulence, the memory-loss by turbulent dis- 
persions is essential. We find a scattering by the small- 
time dispersion and a diffusion by the large-time disper- 
sion in the process of memory-loss and for the derivation 
of the eddy diffusivity. Our formula (32) of eddy dif- 

fusivity can be degenerated to the Heisenberg form (30) 

2 0 

when the scattering prevails , and to the form (31) when 

the diffusion is dominant. The latter form has been used 

as an empirical formula in modeling the atmospheric tur- 

2 1 

bulence by Gisina ; it has also been derived analytically 
2 2 

by Tchen from a cascade theory. The transfer of ens- 

trophy is found to cascade down the spectrum toward the 

2/3 - 3 

large wavenumbers. The spectral law F (k) = Ce k for 
the velocity fluctuations and the numerical coefficient 
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0=2.59, as derived in (51), agrees with the results from 

1 4 

the triad interactions by Kraichnan. The nonstationarity 
of will cause a logarithmic factor due to the inter- 
mittency of turbulence, and will be considered on a sep- 
arate occasion. 
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APPENDIX. SCALING THEORY OF TRANSITION 


We distinguish between a direct transition and a 

m 

retrograde transition of paths Mt) and £(-t) in the posi- 

tive time t and the negative time -t, respectively. The 

probabilities of transition have been investigated by 
2 3 

Tchen , and were found to be governed by the Fokker- 
Planck equations. The derivations were based on a pheno- 
menological ground. Here we shall re-examine the retro- 
grade transition by a scaling procedure and on a basis 
that is consistent with the specified micro-dynamical state 
of turbulence, i.e. the orbit equation: 

-- ^ (~t) . (A1 ) 



see ( 1 4b) . 


The micro-dynamical state of the orbit can be de- 
scribed by the master equation 


(\ + Ljl ) ° , 


(A2) 


upon the introduction of the 6-function 

?(-V,£) = SfH-k-v)] (A3) 


and the differential operator 


L = - tu(-z) ■ • (A4) 
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We may consider the independent variable l to be a random 
variable, so that P becomes a distribution function which 
satisfies the normalization condition 


J 4 *- 1 


(A5a) 


and gives the moment 


hi i p = Xfa) 


(A5b) 


It is easy to verify that the moment of (A2) reproduces 
the orbit equation (A1 ) identically. 

Now we apply the scaling operators A q and A' to (A2) 
and find the following coupled system of equations: 


(\ + \ L i ) K = A, *'(-*) ■ 

- ~&J L 


(A6a) 


(A6b) 


We integrate (A6b) to get 


P'= jlx' . 


21 


and 


(A7) 


A. *(-*)■ 




~>JL 7J2 J ->JL 


(A8) 
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where 


" J 4t' /^(-zJU (-\ u'fz -T')y (A9) 

~ 0 -v /w 

is a diffusivity^ and is a propagator or evolution 
operator in the retrograde transition. The propagator 
U„ may be considered as the inversion of the differential 
operator L^. 

When we substitute (A8) into (A6a), we derive the 
equation of the retrograde transition of the macro-group 




+ 2 _. K' 2A 


(A10) 


in the form of an integro-dif ferential equation. 

In an analogous manner, we derive the equation of 
transition 




+■ 


2 _. 


Here the governing diffusivity is 


(A1 1 ) 



Z 


(A12) 


as caused by the velocity fluctuations except u' . 
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We make a Fourier transformation of (A10) and (All) 


by the formulas 




= A. 


- ck'-Ai-t) 


(•) 


{txf H- T > «) ' A 


-T - t>- J C-«; 


(A13 ) 


where « is an independent variable in the Fourier space, 

•v 

and d is the number of dimensions. We get 


= p y-rj - /r V •• /c; P<f-r, ^ 

^ _l 

«4 - W-' 


( A1 4a) 


Pfa*)- (A1 4b) 


The solutions are: 


4 v -tK»-Xpd , 

Qt-nrJ ” v ° : ty) 


» X.(-x,k) Mjp (- yy :jd r‘ Kj f ‘ ; J 

= */>(■ £ ™ 

= Z('~c,k) w) 


(A15a) 


(Al5b) 


* ✓ 


with the following components of the orbit function: 
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= tjuji (- tV UV ) 

#(-*»)• s') 


(A1 6a) 


(A16b) 




(A16c) 


Note that K. is the diffusivity in a matured (i.e. 
z' 

large-time x) dispersion, and is the diffusivity in 

an unmatured (i.e. small-time x) dispersion. In the 
latter case the orbit belongs to a free flight, so that 


S’M} - / (We/ Ay> «ai 7 

is locally stationary, and the trace in isotropic turbu- 
lence becomes 


% I = ti /«(*' y &\k') * > < r*')y 

/v-v ft* V ^ V * 

= tv ft (u^fa) J ; (A1 8 ) 


the integration over a region of radius k gives 



^ 43ZL hO 3 
J 


in 2 dimensions 
in 3 dimensions 
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Here x is the factor of Fourier truncation. In terms of 
the spectral function F u (k), such that 



in 2 dimensions (A1 9a) 


4-tc f V* )( TJx) 

^ ^ -V v* J 

we find 


in 3 dimensions (A1 9b) 



in d dimensions (A20) 


In conclusion, with the use of (A20) and the introduc- 
tion of the frequencies m^°^ and such that 


^ )Z s j- IV /C * (uf° J <J 0) ) _ £~ Z y 3 p (kO (A2 1 a) 

- V v .v ^ U. S 


and 


m' = w: K'= k’K' 

K ~ - r 


(A21b) 


we find the orbit components to be 
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(A22a) 


tf) - <U^> (- uP )Z "C z J 

'L'C-l, **) - **f> C“ ^ ' 


Now the solutions (A15) can be applied to 


and, subsequently 


A A a -t. 


in-iC-V 


by the use of the formulas (A13). We obtain 




-t ~ ’ v ° 


A * 



and hence 



4(-W ^ OJ (~ v ’ hr ) *) 


(A22b) 


calculate 


(A23a) 


(A23b) 


(A24) 
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CHAPTER D: QUARTER FOUR 
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ABSTRACT 

For the description of the microdynamical state of turbulence, a 
Liouville equation is taken that is equivalent to the basic hydrodynamical 
system of equations. The equation has the advantage of being homogeneous 
and contains less nonlinear terms. Our main objective is the derivation of 
the kinetic equation of turbulence which has a memory in the turbulent 
collision integral. We consider the basic pair-interaction, and the 
interaction between a fluctuation and the organized cluster of other 
fluctuations in the collection system, called the multiple interaction. 

By a group-scaling procedure, a fluctuation is decomposed into three groups 
to represent the three coupled transport processes of evolution, transport 
coefficient, and relaxation. By exploiting the property of quasi-sta tionari ty 
at the different levels of degradation of coherence of the groups, we 
develop a transport theory with the closure by the memory-loss. The kinetic 
equation of the scaled singlet distribution is capable of investigating 
the spectrum of turbulence without the need of the knowledge of the pair 
distribution. 

The exact propagator describes the detailed trajectory in the phase 
space, and is fundamental to the Lagrangian-Eulerian transformation. We 
calculate the propagator and its scaled groups by means of a probability of 
retrograde transition. Thus our derivation of the kinetic equation of the 
distribution involves a parallel development of the kinetic equations of 
the propagator and the transition probability. In this way, we can avoid the 
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assumptions of independence and normality. 

Our result shows that the multiple interaction contributes to a 
shielding and an enhancement of the collision in weak turbulence and strong 
turbulence, respectively. The weak turbulence is dominated by the wave 
resonance, and the strong turbulence is dominated by the diffusion. 
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I. INTRODUCTION 


The microdynamical state of turbulence can be described by the equation 
of evolution of the distribution function f(t,x,v) in the form: 


s\ , A 


(\ f L ) f- - V't . 


(la) 


The differential operator 

■ V -h E-'b c> JS 


L ■ v 


(lb) 


contains a fluctuating self-consistent field E(t,x) which perturbs 


the-, ballistic orbit that is represented by c>. + v'- V in the phase space 
of position x and velocity v • 

A 

In plasma turbulence, (la) is called the Vlasov equation, when E is the 
electrostatic field multiplied by the ratio of electric charge to mass* 

✓A 

The equation can also be used to treat the fluid turbulence, when E 
represents the pressure-gradient or other additional hydro dynamical forces* 

We decompose the. total distribution function 


/N — ~ 

f = f + f 


( 2 ) 


into an average distribution 

f s O > 

by means of the averaging operator A , and a fluctuation 

-A-/ 

f = Af. 

by means of the fluctuation operator 

A = 1 - A . 


(3a) 


(3b) 
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The difference between A and )> lies in that the operator A , 

like any other operators, applies to all the functions which may follow, 
while ( y is an ensemble average of the function or functions limited 
by the angular brackets. 

Thus by applying X to (la), i.e. by premultiplying (la) by A , 
we obtain the equation of evolution of f , as follows: 

(\+L)J~C. W 

The collective collision 

C = - A U , (5) 


can be written in the form 


( 6 ) 


The collision operator 

( 2 ) = 3 •<£> ' f ( 

has a diffusivity ( which may serve as an operator to keep a memory. 

On the other hand, by applying A to (la), we obtain the equation of 
evolution of f . This equation may take the following alternative forms: 


( + A L)f - - L f 

. _ — . ^ ^ 

( ^ + L)f = - L f + C 

O + L) f = - L f - C , 

t 

where C is defined by (5) , and 

/V . — / "V- — f /-%✓ ■'V' 

C = - ALf = - Lf-C . 


(8a) 

(8b) 

(8c) 


( 9 ) 
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Our problem of the statistical theory of turbulence is to find the 
kinetic equation (4), with the turbulent collision (5) to be derived on the 
basis of one of the three equations (8a) - (8c) of fluctuations. The three 
equations differ by their governing differentials 


+ AL 


+ l 


A 

+L 


( 10 ) 


so that their integrations will have different evolution operators, as follows: 


A , 







(ii) 


* 1-3 

The operator y\ was introduced by Weinstock, ^free * s t * ie free 

flight or ballistic propagator, describing a trajectory that is not 

■A A 

perturbed by the E - field fluctuations, and U is the exact 
propagator. 

The differential equation (8a) has the simplest source, but its 
operator is complicate. On the other hand, the integration of the 

differential equation (8b) can be achieved by the simplest propagator 
Uf^ee f ree flight, but the source is complicated by the nonlinear 

fluctuation C , as this will generate a hierarchy and requires a 

4 — 

closure. In the following we choose (8c), because the average source C 

Ay i A 

is easier than the random source C , and the propagator U is simpler 
than the operator /\ , and is governed by the well determined equation 

of the exact trajectory in the form: 


[Z_ + L(t) ] U(t.t') = 0 

or 

L(t')] U(t,t') = 0 


(12 a) 


(12b) 
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(13a) 


~ - 5 

The fluctuating source -L f gives a pair collision: 


C ■ <A>{f } ' 


that is of the diffusive type with an operator 


<A> = 


related to 


the diffusivity (]£)'} 


(13 b) 

The two diffusivities are defined by: 


(i>) = f dZ (E(t,-x).A (t,t-T)E(t-T )^> 

- J o 

(f)) = J dx <Je( 1 :,x) U(t,t-T )E(t--£)^) . 

' J o 


(14a) 

(14b) 


Different levels of approximations have been made in the literature 
on the theories of turbulence in fluids and plasmas- As the lowest 
order, the quasi-linear theory assumed a propagator of free flight 
and obtained a resonance denominator O - k*v, where cj and k are the 
frequency and the wavenumber of a wave 


- i(0 - k-v)T 


(15) 


interacting with the particle of velocity v- This resonance was 
fundamental to the Landau damping in the quasilinear theory- The next 
order approximation assumed an expansion around U in the renormalization 
theory. Interesting advances in theory of turbulence have made use of 
this assumption. ^ 
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With the group scaling, the derivation of the spectrum of turbulence 
can be made on the basis of the kinetic equation of the singlet distributiom 

g 

without the pair distribution. We represent the three transport processes, 
of evolution, transport coefficient, and relaxation, by three scaled 
Croups and analyse their interactions. The closure is obtained by the loss 
of memory in the processes and by exploiting the property of equasi- 
stationarity in the degrading coherence of groups. 

The paper is organized as follows. After a comparison between the 
hydrodynamical and the kinetic descriptions of the microdynamical state 
of turbulence (Sec. II), we introduce a group-scaling procedure for 
representing the three transport processes by the three groups (Sec. III). 

We find that the kinetic equation of the macro-distribution has a turbulent 
collision integral with a memory. We investigate the pair interaction and 
the multiple interaction between a micro-fluctuation and the cluster of 
other micro-fluctuations that form a macro-group. Various forms of memory may 
develop: 

(i) First, the memory .appears in the collision that is controlled by 
the diffusivity acting as an integral operator. Then, since the collective 
collision includes the pair collision and the multiple collision in the 
form of an integral equation, a memory ensues. 

(ii) The Lagrangian representation of a function is made by means of 
a propagator that is governed, by a kinetic equation with a memory. 

Obviously, the memory group (i) which deals with the evolution of 
the distribution function at the time instant t has a longer duration than 
the memory group (ii) which deals with the evolution of the propagator 
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for the shorter time-interval t-t 1 . We shall neglect the memory of the 
group (ii), by justifying from the property of quasi-stationarity 
between the groups (Sec, IV). We shall show that the propagator is related; 
to the probability of retrograde transition, and develop two parallel 
kinetic theories on the basis of the detailed dynamics of the trajectory 
(Sec. V). Finally, the collective collision is investigated, by means of a 
closure based upon the memory-loss (Sec. VI) . A memory function appears and 
is controlled by the resonance function in weak turbulence and by the 
diffusion in strong turbulence, yielding two opposite outcomes: a 
shielded collision in weak turbulence, and an enhanced collision in 
strong turbulence (Sec. VII). 


II. MICRODYNAMICAL STATE OF TURBULENCE ri 

The microdynamical state of turbulence can be described in two ways. 
First we can use the Navier-Stokes equations for the incompressible fluid: 


( ^ + Cl-v - V V 7 ) u. **■ E 


y ■ u, ~ O 


Here u is the fluid velocity, 0 is the kinematic viscosity, and 


(16a) 

(16b) 


E - 7 P 

r ~ 


(17) 


is the pressure gradient, and may even include other forces, such as 
buoyancy and random sources. 


D-8 


Alternatively, we can consider the microkinetic equation (la) , 
rewritten in the form of the Liouville equation: 


[<^+L(t)J f(t,x,v) = 0 . (18) 

A 

The evolution of the exact distribution function f(t,x,v) is controlled by 
the differential operator 

L(t) = v-7- V 7+ E(t,x). e> . (19) 

(V v« 


The microkinetic equation has been used as a basis of deriving the 

EBGK equations by writing f(t,x,v) in the form of a summation of 

8 

£ - functions for N particles in the phase space. Presently, in 

✓A 

order to secure the equivalence of the two descriptions, we write f 
in the form of a single f -function: 


f(t,x,v) = $[v - u(t,x)J . (20) 

The density in the incompressible fluid is taken as unity here without 
loss of generality. The compressible fluid should have the density 
(t,x) as a factor before the ^-function. 

By taking the moments, it can be demonstrated that the 
microkinetic description is transformed into the fluid description as 
represented by the Navier-Stokes equations. 

a 

By the condition of incompressibility, we can write E in the following 
two forms: 

S * - J_ v 

to * J V |x-X'/ 


v'y' * &(%*) u(t,*j Va*/ X > ( 20 ) 

V V y V V ^ V * 
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( 21 ) 


A 


E = 




One is nonlinear and the other is linear, so that the differential 
operator becomes 


L 





with 


(22a) 


l 



(22b) 


The Navier-Stokes equations (16) and the Liouville equation (18) 
may be called the primitive equations in the physical and phase spaces, 
respectively. 

The kinetic approach has several advantages. First, it transforms 
the inhomogeneous Navier-Stokes equation into the homogeneous equation 
(18) and rendered all nonlinear terms, as arising from the velocity 

A 

function u(t,x) , into linear ones with the independent variable v , as 

v « «V 

A 

in (18) and (21) . The only nonlinear term kept is connected with E in 
(18) for describing the mode-couplings. It can ultimately be treated by 
the linear equation of state (21) . 

Secondly, the kinetic approach treats both the wave-wave interaction and 
the wave-particle interaction. The latter is the mechanism of the linear 
and nonlinear Landau dampings or amplifications. This interaction is not 
explicit in the hydrodynamical approach. 
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Finally, should one propose a general kinetic theory that combines both 
the molecular motion and the turbulent motion, the raicrokinetic equation (18) 

«A 

can again serve as a basis, provided E includes both motions above, without 
a predetermined viscosity. Our kinetic theory that is based upon (18) with 
a differential operator (19) is concerned with the turbulent motion in a 
viscous fluid medium. However, for the investigation of the turbulent collisions 
by eddies of size larger than the viscous cutoff, the effect of V can be 
neglected. 


III. GROUP-SCALING PROCEDURE 

The decomposition of a fluctuating function into an average and a 
fluctuation is the usual practice in statistical methods. We denote the global 
average and the fluctuation by 

A and A = 1 - A , (23) 


respectively . 

Not all scales of fluctuations perform the same role in the three 
processes mentioned earlier: evolution, transport coefficient, and relaxation. 
For describing these processes and their couplings analytically, we 
represent them by three scales, using the three operators 

A (o) , A’ , A" , (24a) 

to form the three groups, or scaled fluctuations: 

E° = A° E , E' = A' E , E" = A" E . (24b) 

They fluctuate and have their durations of correlation 
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T c > T 



(25) 


in the increasing degrees of incoherence. This is the property of 
degradation of group coherence, instituting a quasi-stationrity of one 
group with respect to the other. The fields (24b) will also be called 
macro-, micro- and submicro-f ields, respectively. 

The two-scale averaging procedure is obtained by further scaling A 

into 

A = A° + A' , (26) 


by using A^°^ to find 

A° A = 0, A“A“ = A w , A~A' = 0. 

Subsequently, the three-scale averaging procedure is obtained by another 
scaling of A T into 

, _ .(i) 


A' = A + A" , 


(27) 


( 1 ) 


and using A to find 

A (1) A = 0, A^A° = 0 , A^ 2) A (1) = A (1) , A (1) A" = 0. 


The averaging procedures of many scales may be denoted by 

A =< > , A o > o , A 1 ^< / > 1 , (28a) 

where A is the operator of the global average, A q and A^ are the 
operators of the accumulated averages : 

A q = A + A° , A 1 = A + A° + A (1 ^ . (28b) 
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Although the groups A° and A 1 may overlap in their wavenumbers 
instantaneously, their statistical functions, e.g. their spectral contents 
must lie in adjacent domains of wavenumbers, i.e. 

(0, k) , (k,0O) . (29) 

We conclude that the group-scaling procedure is a scaling into the 
three transport processes mentioned above. It can distinguish between the 
two diffusivities: the asymptotic diffusivity ('D') and the non-asymptotic 
dif fusivity ( } , a distinction not present in the one-scale averaging 
procedure. 

Recall that the spectral function, being the Fourier transform of the 
correlation function of velocities at two points, would, in the one-scale 
averaging, require a pair distribution function and its coupling with the 
singlet distribution function. In the present procedure of group-scaling, 
the singlet distribution function f° suffices. Indeed, f° gives u° and 
((u°) 2 > , and thus derives the spectrum of u- fluctuations by a differentiation 

with respect to k 

The degradation of group coherence (25) indicates a consecutive 
quasi-stationarity among the groups, and is an important property for 
classifying the interactions and placing the memory and the memory-loss. 
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IV. TURBULENT COLLISIONS AND MEMORIES 


A. Collision in the Kinetic Equation of f Q 

By scaling the Liouville equation (18) by means of A q and A f , we get 
the equations of evolution of the macro-distribution 

(\+A'L,){=Q oo) 

and the micro-distribution in the form: 

- (31a) 


or, equivalently. 


0 t f ■ < 31b > 

Here 

C « - A L ! f T (32) 

o o N ' 

is the collective collision. The scheme falls into the framework that was 
described by (8c) and (8a) . 

The differential equation (31a) can be integrated to give 

f ' = -A' Jut' Xj(t,t') L'ftJ fjt') 

° f- , (33) 

-A’ [%’ V(t,tj CJf ) . 

^0 
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The lower limit of integration is set to zero, and the initial value is 
ignored, since 


A q L’(t) f'(0) = 0 


( 34 ) 


Here and in the following, the dependence on x or 
The propagator 


u(t.t') 



(x,v) is understood. 


satisfies the differential equations (12) . 

It is to be remarked that the differential equation (31a) for the 
evolution of f T has two sources written in the right hand side: The micro- 
source -L f f o represents a coupling between the micro-operator L f and 
the macro-distribution f Q and gives a fluctuating contribution to the 
solution, as seen from the first term in the right hand side of (33) . The 
macro-source C q represents the collisional cluster in the medium in which f 
evolues. Although is quasi-stationary, it can give a random contribution 

through the operation A 1 !] as seen from the second term in the right hand 
side. This term represents the effect of the quasi-stationary riding on 

the random trajectory. 

In order to determine the collective collision, as defined by (32) , we 
multiply (33) by L f (t) and average, obtaining: 


C (t) = A dt 1 L ' ( t) U(t, t 1 ) L’(t') f (t f ) 
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+ A dt' L' (t) U' (t,t') C (t’) 
o / o 


( 36 ) 


This expression consists of the pair collision 


<C°(t) = A / dt’ L* (t) U(t,t') L'(t') f (f) , 

y 0 


(37) 


and the multiple collision in the form: 


A o / dt* L'(t) U'(t,t') C 0 (t') = 


- A o / dt' H o (t,f) C 0 (f) 


= - A H * C , 
o o o’ 


(38) 


Here H can be identified as the collective collision in the kinetic 
o 

equation of U q , i.e. 

£ ^ + L o (t) ] U o (t ’ t?) = _A o L,(t:;> " H 0 (t»t’) , (39) 

as obtained upon scaling (12a) by A q . Note that (12b) will not lead to 
the desired collision and is therefore not relevant here. 

As an option to (8a)-, we have the equation of evolution of the 
micro-distribution in the form 



(40) 


and obtain 


f' = - A* / dt'/l(t,t') L'(t') f o (t») 


(41) 


D-16 



by an integration and the use of the operator /\_ . follows: 


C = - A L'f ’ 

o o 


dt' L*(t)A(t,t') L'(t') f 0 (t') 


A o A- {*„} 


( 42 ) 


Here A q /1 denotes the collision operator 


A q 4' = A q /dt’ L*(t) A(t.t’) L’(t') 


= c> • A q D 1 • 'b , 


(43) 


with a diffusivity 


AD' = A / dt' E'(t) A(t,t') E'(t') . 

o o * 


O -a 


(44) 


By replacing by U in (42) - (44) , we obtain the pair collision 

C “ A o A ' { f o } • <45 > 


with a collision operator 


A q A = A q / dt' L'(t) U(t,t') L'(t') 


= ■ A. IP ' 4 ) , 


(46) 
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and a diffusivity 


A D' = A f dt’ E'(t) U(t,t’) E'(t') 


(47) 


By collecting the results (36) - (38) , we obtain the following 

relation between the two collisions C and £T in the form: 

o o 


C « f ~ A H *C 

o % 000" 


(48) 


wi th 


A H *- C ■ « A 
000 o 


dt’ H o (t,t T ) C Q (t f ) 


(49) 


B. Collision in the Kinetic Equation of U q 

The distribution function f and the propagator U have their evolutions 
of some resemblance. By repeating the calculations that have been made in 
Subsec. IVA, we obtain the collisions 


H 

0 

- 

O 

/ 


(50a) 

*0 

= A A 1 
0 

[" 0 } 

(50b) 

U , as related by 
0 

the formula: 


H 

0 

0 

a 

11 

- II * H 
0 0 

(51) 


Here 11 and ll| are the collective collision and the pair collision 
0 0 

in the equation of evolution of U q , i.e. (39). The results (50a), (50b), 
and (51) are analogous to (42), (45) and (48), respectively. 
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C. Memories 


We summarise the results obtained above as follows: 

(i) Primary memory in C q 

The kinetic equation of f , in the form (30) and with the collision 


(42), carries a primary memory by the collision operator 
be written into two groups as follows: 

j, which can 


(52a) 

The second group forms a correlation \ on its own, but does 

not have sufficient time to approach its equilibrium, by (25) , and therefore 
will not contribute to a diffusive collision. Hence we can write 

- Wif } 

(52b) 

as an approximation, reducing (42) - (49) into: 


c.-<*‘W ' i} 

H -<A>-0 , H -<A'> o 

(53a) 

(53b) 

C. = C o - H * Q 

(53c) 

f - ft - h *h , 

(53d) 


with 

(O- t ■(¥>■* , <A'> = v 

(54) 

6 ') -.j%'(EXt)A(t l t'jt(t')y (ip') 'f4t '(e ‘(t) O(Kt’) E(t’/) 

3 4 ) v " ~ h 
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Hence the evolution 


( V A oV f o (t>3 !>! ) = < A ^o (t,) } 

/• t 

- J (i '(t) A (t t-l) L '(t'tpfz 0 (t-TT >J. , (55) 

with t-t 1 =T’ , carries a primary memory in from the fluctuations E 1 

of scale 77 > such that 
c 

t ? \ > T , (56) 


by (25). 

(ii) Secondary memory in H q 

On the other hand, if we write U(t,t f ) = U(T ) , with t-t T =17, 
the kinetic equation 

(\+aL 0 )u o (t, t-t ) = H 0 (t, t-vj 

= /V (lYv)A(tA-^)LYt-z')) VJv-t'J (5?) 

J o 

has a collision integral governed by the collision operator 

f '(L'(r )A (t / ~c-~c , )U(x~t , )\ 
v 0 ' 7 

V 

+ \ Z Jix\L[t) C(t-xp 

u 

decomposed into two parts. With the condition of quasi-stationarity 
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as a continuation of (56) from (25) , the second integral has a matured 
diffusivity, where we can put X-^ce> and obtain while the 

first integral has an unmatured diffusivity. Since T X » the secondary 
memory carried by the collision integral H q in the short time span T 
can be assumed negligible as compared with the primary memory. 

D. Resolution of the Integral Equations 

From the definition 



-e 




HM, 


the Fourier transform of the convolution H *- C is 

o 


JL 

7C 




Hfc') C 0 (t-~c'J 



r oo 

7 4V -e 


_Ccot: '' -lu)t 



(59) 


and is written as 


H #C 0 



(60) 


Here a change of variables ~C~~C has been made. The symbol r=r denotes 

a Fourier transformation. Hence we transform (53c) into: 

=Q(uj) - 7 vHfaJC/aJ. 


(61) 



■■ HI 


V. KINETIC EQUATION OF THE TRANSITION PROBABILITY 


A. Path Dynamics 


The trajectory with which the propagator evolues is described by 
the following differential equation of the path dynamics: 



or, identically, by the system 


d x (t'J 




with the conditions 

x(t) = X, 


v(t) - V , 


which specify that at the instant of time t , 'the said trajectory 
passes by the phase point (x,v) . 

By integrating (63) and writing t 1 = t 9 we find the dynamical 

variables 

*(t- -D) = A + y \/(t-T.) = V -h V('T) 


and the displacements 



(62) 


(63) 


(64) 


(65) 


(66a) 


(66b) 
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Now the path dynamics (63) that is driven by E has the displacements: 


d r 

d v( t-z) 
4x 


rr - V 




= - £(t-v) 

V* ' 


(67a) 


(67b) 


B, Retrograde Transition Probability 

The Lagrangian correspondent of an Eulerian function G(t,x) can 
be written in the form 


<M 

G jTt-T, x+ ^(-T) I t,x J , withr>0 , 


(68a) 


and specifies the value of the function as observed by a fluid particle 
coming from the point x(t -T ) at the time t- ~C > along the 
trajectory which at the time t passes the point x . The Lagrangian 
function can be written more conveniently as 

U(t,t- V ) G( t-T ) (68b) 

a 

by means of the propagator U(t,t-*k ) . It is in reality a function 

of two states: 


t -T » x + Ji (- T ) and t , x 


(69) 


/A 

Since $ (-"C ) is a random function, we can introduce a retrograde transition 

0 

probability, written in the form: 
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111 


p(t-T; xi-l , 


(70a) 


or, briefly in the notation 


(-'C, -A y 


(70b) 


The first form expresses the state t,x as a condition and can be called 
a conditional probability. The second form implies a quasi-stationary 
transition, i.e. p varies more rapidly with — than with t,x . 

-v — w 

Here Jt is an independent variable. 

The concept of the retrograde transition probability has been introduced 

earlier and its equation of evolution was found to be of the Fokker-Planck 
9 

type. More recently, by the use of this transition probability, the 
kinetic equation of turbulence was found . ^ In the present Section, 
we shall devote to the relation between the retrograde transition probability 
and the propagator. 

As the basis of the dynamics of we write the equation of 

the detailed evolution in the form 


+ i(t-T) ) = O , 


(71a) 


satisfying the condition of normalization: 


Jail Jo(-T y .i) - I 


(71b) 
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The differential operator is 


Z(t-v) 





( 72 ) 


In order for (71b) to describe the microdynamical state of the turbulent 
trajectory, we put 


ft (-T, L) = S[i - i(~-cj ] . 


(73) 


The partial differential equation (71a) has its characteristical 
equations coinciding with (67) . Evidently, by substituting (73) and 
integrating (71a) with respect to Jl , we will reproduce the basic 
dynamical equations (67) . 

The Fourier transforms of (71a) and (71b) are: 


(Tc' h ft)] ft (~ T '^ ) = ° 

and 

p(-T,k = °J = (ztvJ ^ , 

with 

2(t-T,k) = -Ck-y(t- v) 


(74a) 


(74b) 


(75a) 


and 


J 4t' k) = -ik-M-tJ 


(75b) 


Here d=3 in the three dimensional space. The integration of (74a) 
with the condition (74b) gives the solution 
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^)(-X,k) = (zit) Jmp [~f c/v t i(-t'kj] 


t -4 


•= (Zic-J -eoyti 


(76a) 


He v 

It consists of a ballistic orbit function e — and a field-dependent 
orbit function, which is 


/v 

T(-T,k) = **/> ik'A(-t) , 


(76b) 


if E = 0 


C. Relation Between the Propagator and the Transition Probability 

The evolution o^- U is governed by a kinetic equation with a 
^ 

collision H . It is to be noted that U(t,t ) , being an evolution 

operator, or propagator, describes the exact trajectory, and therefore 

A 

is a functional of the path during the interval of time from t 

to t-TT in the retrograde transition. Hence the propagator can be written 
as : 

U(t,f') = , (77a) 


or 


U(t.t') 



u (i ) P(-T J), 


t-t' = v , 


(77b) 


by the use of the transition probability p(-T , ■* ) - In an analogous manner, 
we have: 


vat; = J4 vttj p(-*M ■ 


(77c) 
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In this way, the propagator is determined by the probability of transition 
» oe P (— » £ ) % 

D , Lagrangian-Eulerian Transformation 

The Lagrangian function depends on two states, and can be written in 
three forms, by means of a functional as in (68a), a propagator as in 
(68b), and now a transition probability in the form: 

U(t,t-D G(t-T ) p(-T , /) G(t--C , x+/f) . (78a) 

It admits a Fourier form: 

U(t,t-T) G(t-T)'=r (2^c) d p(-T,k) G(t-T,k), (78b) 

by a Fourier transformation with respect to X • A time integration 
yields the convolution: 


U * G 


f 

V c\ 


dz u(t,t-r ) G(t-r ) 


(79a) 


Hence the combined Fourier transformations with respect to both 
X yield •' 

O’ * Cr (z-ji 1 x CJ,k ) (r(ud,i ij 


or 


O * Cr - K [J J-r ] Cr(u,kJ ■ 


t and 


(79b) 

(79c) 
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Use of (76b) has been made. 

We conclude that the transition probability helps in determining 
the propagator and the function to be operated upon, and thereby establishes 
the Lagrangian-Eulerian transformation. 


E. Kinetic Equation 

By scaling (71a) by A and X , we find the system: 


= J(- T ‘V <80a> 

^ -K-tA], (80b) 


with the collision 




(80c) 


The procedure of determining the collision is the same 

«v 

as that yielded the collision H in the kinetic equation of U .We 
deduce 

J - f - H a J , (81) 


analytically, as was with (53c) and (53d) . We have 





(82a) 
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(82b) 


-*• wv ^ ^ 

with 

<k> =/^r - (y(^j 

(I K ) =/*' (VT-W Cf-T, -T-rfo(-T-T'/)> . 

■z 0 »• 


(82c) 


(82d) 


Now by neglecting the secondary memory, we can write the approximations: 


H = HJ , J = J7 (83a) 

(lK(T } - (iK(-c )) • (83b) 

~v- 


In the following, the determination of U and U will be 
performed through 'p and p~ , so that we will have no more opportunity 
of dealing with the kinetic equation (57) of the propagator. 


VI. DETERMINATION OF THE COLLISIONS BY MEANS OF THE TRANSITION PROBABILITY 

The result (79c) has transformed the convolution U+G of the Lagrangian 
form into its Fourier Eulerian correspondent. We shall apply this result 
to transform the convolution H *• C , rewritten as 

H*C o <A)u*C o * 

by definition (53a). We obtain: 
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H * Q % x <£') Hr f(. x k) Q(o>,k) 

J o - 

= - 7 c (l/,) k z M(o)jk ) C a C^,h) 


where 


= - Q,(u,k) > 


xfakj = T <J>') M(u,k) 


is a dimensionless factor of memory, and 

, 2. r® - k ■ v) t 

* / dv * " 

* J o 

«• 1/3 

is called the memory function", defining a life-time M for the 

memory . 

Rigorously speaking, = operates on the ballistic 

orbit function (15) , as well as on all the dif fusivities which are 
embedded in U and C . With the approximation that varies 

O JV 

slowly with v , as is true in strong turbulence, we can evaluate (85b) 
and get 

mm; a J\ T 1 * (u ' k ' vh W^il 


where, by definition (76b), we have 


(84) 

(85a) 


(85b) 


(86a) 
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(86b) 


?(--£, I ) a CkJj-vj) 

■m . iW mi 

Now it leaves us the function p(-T ,k) , and we shall determine it 
by means of the Fourier transform of the kinetic equation of transition 
in the form 

lp(-T,&) = ffakj (87) 


from (80a) , The collision is given by the expression 

f 

j(-r,k) - - k l {kto) , W8> 

as governed by the diffusivity (82c) • The integration with the condition 
(74b) gives the probability 


p(-t; , k) = (27 V ) d exp k 2 J dr 


(89a) 


which we substitute into (86b) to obtain 

jv 


P(--C,k) = exp^-k 2 J dz^K(v'))]. 

J o 

Since ( is caused by the field fluctuations 
E = E° + E* and E' = E^+ E" , 


(89b) 
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we can decompose the diffusivity into components, as follows: 


(k) = ( k °) + (k^>, ( k ') = <K ( ^+ (90) 


Within the available time span T and under the quasi-stationarity 
condition (58), the diffusivity reaches a matured development, 

leaving the non-asymptotic dif fusivities of the order of magnitude: 


(k ( 1 ) (r)> « <k°(x )) . 

The negligible (t ))> may be included under ^K("ir ) ^ for brevity. 

The procedures indicated are to calculate the two dif fusivities from 
the path dynamics (66b), to substitute the results into (89b), and 
subsequently into (85a/ and (85b), to obtain P , and M . The calculations 
are lengthy, but can be simplified if we choose the highest power in T* 
for retaining the most rapid memory-loss. The results are collected 
as follows: 

(a) results relating to 


h *- C e = -xfak ) C D (^,k) 


?hi) * 




1 : = t-t' 


* (io,k) = c u ^ 3 Mfa,k) . <y - 


(91) 


bfak) * n 


with 


*r r e 




. / / ] ,'37 o ¥- f- 


A 
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(b) results relating to {])'/> 


with 


* C = -tifai) C/cOyij 

<x(u,k) = C, . 

Ml (u,i) * fir t* -t 

Jq 


C(u>-k-*h -<4> 3 t 3 


( 92 ) 


By applying these results, we transform (53b) and (53c) into: 


C 0 (oo y k) = c o (u,k) + QiUj 
. /f^4J -a tfakJHfak) 


(93) 


(94a) 


The last equation can be used to derive 


(94b) 


VII. SHIELDING AND ENHANCEMENT OF THE COLLISION 

With the approximation of neglecting the secondary memory, as stated 
in (83), we can reduce (93) into the form 

Cj«,k) = (C (tyk) 4- (to, kj 

rewritten as 

C/to, k) =■ &) € o (u,kJ f (95a) 
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by introducing the coefficient 

(i(w, 1 6 - )) 

- [(!-«/* 

where $(, are the real and imaginary parts of 
'/ *■ 

t X = &! + i <X Z • 


(95b) 


(95c) 


The result (95a) shows that the collective collision depends on the 
factor ^ , which in its turn depends on the memory function jrfj . This 

function has a complicate integral (92) . We shall calculate it 
approximately by means of an interpolation which correctly covers the 
three regions as dominated by UJ , m , and k- v , separately. 

For the sake of simplification of writing, we introduce the following 
frequencies : 


r = 

JPtrv 




(96a) 




CO 


~k-vj 

-V '**' ' 


with the ratios 


5 = C/t ' 1 - t /t 

' Jfin, / D Tt> 




(96b) 


and the numerical coefficients found to be: 


VO&j] 


/ 




(96c) 
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Finally we find the results, as follows: 


<*, = 7' 






- 1 


(i-3?*) L + v6+*ry' ( 97 ) 

= ?*«, • 

We estimate that yj is of the order of unity, and ^ is an index of the 
strength of turbulence. The asymptotic cases reduce the general formula 
into the following: 

(a) For weak turbulence, we have ^ >j> / , entailing we 


find 


so that 


K 


* = (Va) 


(98a) 


(i> - O + X/J 


- / 


(98b) 


causing a shielded collision. 


(b) For strong turbulence, we have 9 entailing Gf <<£. , we 

^ / 


find 


]'< l 




(99a) 


so that 


ft = 0- x,) 


-2 


(99b) 


causing an enhanced collision. 
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can be evaluated for 


The enhanced collision by the factor ^3 

plasma turbulence and fluid turbulence. As an illustration , we consider 

the plasma turbulence in a strong magnetic field. Here the electric field 

fluctuation excites the velocity fluctuation. We found that both spectra 
—3 i ' 

possess the k power law, that ^ is a constant quantity, and that 
the Larmor radius is a spectral cutoff. Since the E-fluctuation drops 
rapidly, the macro-electric energy ^(E°)^ will cover the main body of the 
energy-containing portion of the spectrum to become approximately 
independent of k . These parameters determine a turbulent Reynolds number 


R e - 1/2 


( 100 ) 


By definition (96a) , we have 




obtaining 




and 


|9 * I + *(Vif [i+ n*(lA )'] 




(101a) 


For high intensity turbulence , is large, and we reduce (101a) 


into: 



+ 



with 


(101b) 
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The collective collision is increased above the pair collision by 

-3/2 

an amount proportional to k . This means that the big eddies in a bath 
of smaller ones will be more apt to have an enhanced collective collision. 

An analogous increase of the effective diffusivity of a suspension of small 
particles in Brownian movements is well known. The application to fluid 
turbulence will be given at a later opportunity. 

VIII. SUMMARY AND DISCUSSIONS 

The point of departure of the statistical theory of turbulence is the 

description of the microdynamical state of turbulence, either by a system of 

hydrodynamical equations in the t,x^ space, e.g. the Navier-Stokes 

equations (16), or by a consistent microkinetic equation in the t,x^,v^ space, 

e.g. the Liouville equation in the form (18). Here by enlarging the 

dimensionality from the passage of the t,x^ space to the t,x^,v i space, we 

have eliminated the nonlinearity as connected with the velocity field, 

but have kept the nonlinearity as connected with the E^-field from (21) . 

If the Navier-Stokes equations are further decomposed into Fourier series, 

the multi-dimensional variable u , with o( g (i,k) , has the three 

<K 

directions i“l,2,3 and all the Fourier modes k running f rom -co toQ3 . 

The distribution 

( 102 ) 

transforms the Navier-Stokes equation into a corresponding Liouville 


f(t -\ > - u « (t) ] 
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equation of the forms 


O t + L)f 

where 

L — 2 $ 

is the differential operator, and V , A ... are the coefficients. 

<K <*(11 

In this kinetic representation, the enlarged dimensionality by the 
Fourier transform renders L deterministic, so that the Liouville 
equation becomes linear. Since the linear equation does not distinguish 
between the average and the fluctuations, it does not lend to clearly 
express the coupling between the modes in a functional form. An external 
fluctuation response assumption becomes necessary for establishing the 
physical and functional structure of L . For example, when L and f 
are decomposed into two components as 

L = (L) o +£ L, f= (f) o +Sf, (105) 

(103) can be transformed into 

[^_+ (L) 0 J (f) c * -<^L If) . (106) 

An iteration along the direct interaction approximation was to assume a 
dissipative (L) q of the Fokker-Planck form. 

By contrast, the nonlinear Liouville quation (la) retains the mode-couplings 
which are described by the equations of evolution (8a)-(8c). The form 
(8a) gives the collision in the most direct manner but requires the knowledge 


0 , 


(103) 


+ Z — /i.yV'fl 


(104) 
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of the operator /[ whose physical structure and relation to U can only 
be obtained through other two equations, (8b) or (8c). In the form (8b), the 
memory disappears in the ballistic orbit and reappears as a high order 
correlation 

<L(t) u free (t,t-r) L(t-r ) f(t-T)> , ( 107 ) 

causing a hierarchy. Therefore we have chosen the form (8c) which does not 
present this difficulty. On the contrary it establishes the necessary 
relation between A and U / thanks to the transition probability. 

For the sake of simplification of discussion and abbreviation, we 
use the one-scale average by writing: 

f(t,x,v) a f(l), f(t,x,v) - f(l) , (108a) 

and 

L(t,x,v) f(t,x,v) = f(l') f(l) , (108b) 

with 

f(l') -£•{?}' , (108c) 

from (22b). The time derivatives along the trajectories will be written as: 

d t «> t + L“ , "d t ■*£> + L . (108d) 

The average distribution and propagator, and their fluctuations, 
are governed by the following equations of evolution: 

\ £ (1) - C(l) (109a) 
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(109b) 


d fc 7 ( 1 ) = - 7 ( 1 ') f(l) - C(l) 

and 

d t u(i) = 5(1) 

d t U(l) = - f(l') U(l) - H(l), 

with the collisions: 

C(l) = -<f(l') f(l)> 

5(1) - -^f(l') U(l)^> . 


(110a) 

(110b) 

(Ilia) 

(Hlb) 


The evolution of f(l) is governed by a variable source coupling 

/v __ __ 

f and f and a deterministic source -C which represents the cluster of 
fluctuations in the diffusive medium. Upon integrating (10%) , we get 

7(1) » - Xu * [7(1') 7(1) + C(l)] , (112) 

and subsequently upon multiplying by f(l') and averaging, we find 
the collision, as follows: 

5(1) -<^7(1') AU# fU*)^ 7(1) + {7(1’) u)^c“(l) , (113a) 


or 


C - - H jr C 


(113b) 


Here 


€ * 


( 114 ) 


is the pair collision, having an operator 


<A> 


fd») au * f(i'r> , 


(115) 


and -HfrC is the convolution of two collisions: 


C \ and H -<2> ; U . 


(116) 


In terms of the collision operator 

/~ 


(& y = (la*) ZA* f(i f )} , 

we can rewrite (113b) as 


(117) 


-P^/ -( 2 )W 2 M F l 


(118) 


We conclude that the variable source and the deterministic sources in (109b) 

yield a collective collision with two components: one is proportional to 

* ant * t * ie other is proportional to ^ , as re l ated to the 

pair collision and multiple collision, respectively. The latter contributes 

to a shielding of the collision in a weak turbulence as governed by 

the resonance of the wave exp £-i(««J - k- v)T J , or to an enhancement 

of the collision in a strong turbulence as governed by the diffusion 

parameters of relaxation frequencies CJ ' and m° in the memory function. 

J* 

Such opposite effects can find their analogies in the Balescu-Lenard equations 

13 14 

in quiescent plasmas ’ and in the effective diffusivity of a syspension of 

particles undergoing Brownian movements, respectively. More recently, 

a reversal of the effect of multipe collision in Brownian motion has 

also been found, dependent on the relative scales of the particles and 

16-18 

the Brownian motion. 
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Ill 


We can give the following interpretation, A cluster of large size a may 
move with a high velocity and be scattered by the small scale turbulence, 
i.e. ka^ 1. Since a part of the eddies will be displaced or shielded from 
participation in the collection interaction, the effective diffusivity 
is reduced. On the other hand, if ka<^ 1 , the cluster will participate 
together with the background eddies to contribute an enhanced collision. 

By the generalization to the two -scale averages, we obtain the kinetic 
equation of the distribution function f with the collective collision 
C Q , as found in (30) and (48). This can derive the transfer function for 
the energy-cascade and determine the spectrum of turbulence without the 
intermediary of the pair distribution function ^f(l’) 
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